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1.  INTRODUCTION 


Modern  electronic  systems  must  frequently  operate  in  hostile  electro¬ 
magnetic  environments  dominated  by  unwanted  signals.  The  unwanted  signals 
are  referred  to  as  interference  whether  or  not  they  cause  unacceptable 
system  performance.  Undesirable  electromagnetic  interactions  can  occur 
between  different  1)  systems,  2)  subsystems  within  a  system,  3)  equipments 
within  a  subsystem,  and  4)  components  within  an  equipment.  Electromagnetic 
compatibility  (EMC)  exists  when,  even  in  the  presence  of  the  unwanted  signals, 
the  equipments,  subsystems,  and  systems  perform  together  in  an  acceptable 
manner.  In  other  words,  if  unintentional  electromagnetic  radiations  and/or 
responses  do  not  sufficiently  degrade  system  performance,  EMC  is  said  to 
apply.  On  the  other  hand,  electromagnetic  interference  (EMI)  results  when 
acceptable  performance  is  prevented  because  of  the  interfering  signals. 

The  term,  interference,  is  seen  to  have  a  dual  meaning  because  it  refers  to 
both  the  unwanted  signals  and  their  effect. 

In  order  to  determine  acceptable  performance,  a  performance  criterion 
is  needed  which  provides  a  standard,  rule,  test,  or  measure  for  judging 
quality  of  operation.  Depending  upon  the  specific  application  and  the  type 
of  system,  subsystem,  or  equipment  involved,  different  criteria  are  appro¬ 
priate.  Several  examples  are  cited  below.  In  speech  communication  systems, 
where  intelligibility  is  the  important  factor,  two  widely  used  criteria  are 
articulation  score  and  articulation  index.  In  television  and  graphics  dis¬ 
play  systems,  resolution  is  often  chosen  as  the  performance  criterion.  Bit 
error  rate  and  probability  of  error  are  commonly  used  criteria  in  digital 
communication  systems.  Performance  criteria  used  in  analog  data  systems 
are  mean  square  error,  peak  error,  and  average  absolute  error.  With  respect 
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to  radar  systems,  the  probability  of  detection  for  a  specified  false  alarm 
probability  is  a  frequently  used  criterion.  In  every  case,  once  a  perfor¬ 
mance  criterion  has  been  chosen,  acceptable  per f ormance  can  be  specified  by 
selecting  an  allowable  range  of  values  for  the  performance  criterion.  This 
range  is  usually  based  upon  mission  requirements. 

In  the  evaluation  of  system  performance,  the  desired  and  interfering 
signals  are  usually  characterized  in  terms  of  convenient  parameters  such 
as  average  power,  peak  amplitude,  energy,  and  peak  power.  The  EMI  perf orm tru-e 
curve  is  a  plot  of  the  selected  performance  criterion  versus  some  function 
of  the  desired  and  interfering  signal  parameters.  For  example,  in  Figure  1 


Figure  1  EMI  performance  curve  for  a  speech  communication 

system. 
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is  shown  a  typical  EMI  performance  curve  for  a  speech  communication  system. 

For  this  system  the  performance  criterion  selected  is  the  articulation  score, 

AS,  which  gives  the  percentage  of  words  correctly  interpreted  in  a  listener 
panel  test.  The  average  desired  signal  power,  S,  and  the  average  interfering 
power,  I,  are  the  desired  and  interfering  signal  parameters  suitable  for  this 
plot.  The  EMI  performance  curve  consists  of  a  plot  of  the  articulation  score 
versus  the  signal-to-inter f erence  ratio,  S/i.  Use  of  the  signal--to- interference 
ratio  implies  that  performance  is  independent  of  the  absolute  levels  of  S  and 
I.  Rather,  it  is  their  ratio  which  determines  the  amount  of  intelligibility. 

As  would  be  expected,  the  articulation  score  is  seen  to  be  an  increasing 
function  of  the  signal-to-inter ference  ratio  until,  for  large  enough  values 
of  S  relative  to  1,  essentially  perfect  intelligibility  is  achieved. 

Given  the  EMI  performance  curve,  it  is  possible  to  specify  the  range 
of  operation  over  which  acceptable  performance  occurs.  The  performance 
threshold  is  that  value  of  the  performance  criterion  which  demarcates  the 
regions  of  acceptable  and  unacceptable  performance.  These  regions,  of  course, 
are  determined  in  accordance  with  mission  requirements.  By  way  of  example, 
the  performance  threshold  in  Figure  1  has  been  set  at  an  artic.nl  ation  score 
equal  to  60%.  This  implies,  for  the  application  of  interest,  that  acceptable 
performance,  results  when  the  articulation  score  is  greater  than  or  <  qua l  to 
60%  while  scores  loss  than  60%  yield  unacceptable  performance.  To  put  it 
another  way,  F.MC  exists  provided  AS  _>  60%  while  EMI  prevails  when  AS  <  60%. 

Given  the  F.MT  performance  curve  and  numerical  values  for  the  performance 
threshold  and  the  Jesired  signal  parameter,  it  Is  possible  to  specify  the 
level  of  interference  above  which  EMI  exists.  For  example,  with  respect  to 
Figure  1,  S/l  equals  -18  dB  at  the  performance  threshold.  Therefore,  when 
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S  =  -10  dBm,  Interference  with  average  power  greater  than  8  dBm  produces  EMI. 
On  flu!  other  hand,  with  S  =  -10  dBm,  EMC  exists  when  I  £  8  dBm.  By  definition 
the  suscept ihillty  1  eve  1  ,  L,  is  the  level  of  the  interfering  signal  th.it 
results,  for  a  specified  desired  signal,  in  operation  aL  the  performance 
threshold.  Although  the  numerical  value  of  L  depends  upon  the  desired 
signal, the  susceptibility  level  is  primarily  a  characteristic  of  the  system, 
subsystem,  or  equipment  involved.  Returning  to  the  example  of  figure  1, 
assume  S  =  -10  dBm.  The  corresponding  susceptibility  level  is  given  by 
L  =  8  dBm.  However,  an  improved  design  of  the  speech  common ic. at  ion  system 
of  figure  1  could  result  in  the  susceptibility  level  being,  increased  to  a 
value  greater  than  8  dBm  when  S  =  -10  dBm. 

It  is  important  to  recognize  that  EMI  performance  curves  are  a  function 
of  the  type  of  interference  encountered.  For  example,  given  an  KM  communi¬ 
cations  receiver,  separate  EMI  performance  curves  apply  to  the  dixiorer.t  ca  ■ 
of  AM,  FM,  and  pulsed  radar  interference.  Although  a  complete  set  of  KM I 
performance  curves  does  not  exist,  many  curves  have  been  generated  for  sewra 
types  of  communications  receivers  under  a  wide  variety  of  interfering  signal: 

Having  introduced  the  concepts  of  EMI  performance  curve,  performance 
threshold,  and  susceptibility  level,  attention  is  now  focused  on  EMC  modeling 
and  analysis.  Historically,  approaches  to  this  problem  proceeded  Irma  a 
deterministic  point  of  vio.w  in  which  the  signals,  coupling  paths,  and  equip¬ 
ment  characteristics  art.-  assumed  to  be  known .  Uncertainties  were  accounted, 
for  by  utilizing  "worst-case"  models.  However,  many  nr ..  .  rn  elect  r  -  emetic 
environments  and  electronic  equlpme  .  hav,.  !  ecomo  so  complex  tne>  no  longer 
can  be  modeled  cither  realistically  or  efficiently  using  a  deterministic 
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approach.  This  is  particularly  true  of  equipments  containing  microelectronics. 

In  particular,  digital  circuits  are  currently  being  designed  with  the  latest 
integrated  circuit  technologies  where  hundreds  and  even  thousands  of  active 
devices  are  fabricated  on  tiny  silicon  chips.  The  resulting  circuits  are 
so  complicated  they  cannot  be  analyzed,  or  even  simulated  on  a  computer,  in 
order  to  determine  their  EMI  performance. 

In  recent  years  macro  models  have  been  introduced  for  simulating  com¬ 
plicated  circuits.  Ideally,  the  behavior  of  the  macro  model  closely  approx¬ 
imates  that  of  the  original  circuit  even  though  the  macro  model  employs 
significantly  fewer  active  devices.  Limited  success  has  been  achieved  thus 
far.  For  example,  EMI  in  an  operational  amplifier  containing  25  transistors 
has  been  successfully  predicted  by  a  macro  model  employing  only  2  transistors. 

In  addition,  macro  models  are  currently  being  developed  for  predicting  EMI 
in  simple  digital  logic  circuits. 

However,  the  electromagnetic  environment  and  equipment  susceptibility  are. 
In  reality,  random  in  nature.  By  utilizing  a  probabilistic  approach,  statistical 
macro  models  can  be  developed  which  simplify  the  problem  of  characterizing 
complicated  signals  and  circuits.  In  such  an  approach,  detailed  circuit 
models  are  replaced  by  statistical  models  where  probability  density  functions 
are  used  to  evaluate  probabilities  and  statistical  averages  associated  with 
the  response.  Following  the  lead  of  Dr.  Capraro  of  Rome  Air  Development  Center 
[2],  this  report  proposes  and  analyzes  a  probabilistic  EMC  model  which  is 
useful  when  the  deterministic  approach  to  EMC  modeling  and  analysis  is 
inappropriate.  Additional  material  relevant  to  this  work  can  be  found  in 
the  paper  by  Bossart,  Shekleton,  and  Lessard  [3], 
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2.  THE  PROBABILISTIC  MODEL 

For  simplicity,  consider  the  EMC  problem  illustrated  in  Figure  2 
where  both  a  desired  signal,  s(t),  and  an  interfering  signal,  i(t) ,  couple 


Figure  2  A  simple  EMC  problem. 

into  an  equipment.  If  the  desired  signal  is  to  convey  information,  it  must 
be  unpredictable.  In  fact,  according  to  information  theory,  the  information 
content  of  a  signal  increases  with  its  uncertainty.  Unpredictability  in  the 
desired  signal  is  also  introduced  by  the  coupling  path.  Typically,  the  re¬ 
ceived  signal,  s(t),  must  be  distinguished  from  the  transmitted  signal,  d(t)  . 
For  example,  unpredictable  signal  distortion  may  arise  when  the  attenuation 
and  phase  characteristics  of  the  coupling  path  are  unknown.  Also,  the  trans¬ 
mitted  signal  may  undergo  fading,  whereby  the  received  signal  amplitude  is 
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found  to  fluctuate  in  a  random  fashion.  In  addition,  additive  noise  may  be 
introduced  in  the  coupling  path  so  that  the  received  signal  actually  consist 
of  the  sum  of  signal  plus  unknown  noise.  In  the  probabilistic  approach,  un¬ 
certainty  in  the  desired  signal  is  incorporated  by  treating  s(t)  as  a  sample 
function  from  a  random  process. 

The  interfering  signal  may  enter  the  equipment  directly  through  the  input 
port.  It  may  also  penetrate  the  equipment  indirectly  by  coupling  on  to  cables 
such  as  power  lines  or  by  propagating  through  apertures  in  the  case.  As  with 
s(t),  the  interfering  signal  i(t)  is  likely  to  be  unpredictable.  Some  of  the 
uncertainty  may  be  due  to  the  transmitted  signal,  j(t),  also  being  an  infor¬ 
mation  bearing  signal.  Additional  uncertaint  may  be  introduced  by  the 
coupling  path.  Finally,  the  mechanism  by  which  i(tj  enters  the  equipment  may 
not  be  completely  known.  This  is  especially  true  when  the  interfering  signal 
couples  onto  wires  or  cables  and/or  through  apertures.  As  a  result,  i(t)  is 
also  treated  in  the  probabilistic  EMC  model,  as  a  sample,  function  from  a 
random  process. 

The  final  uncertainty  arising  in  the  EMC  problem  is  associated  with  the 
equipment  itself.  It  is  not  possible  to  manufacture  "identical"  resistors, 
diodes,  transistors,  etc.  having  precisely  the  same  parameter  values.  Porasitii 
elements,  whose  values  cannot  be  controlled,  may  play  an  important  role  in 
an  equipment's  behavior.  The  positioning  of  wires  and  cables  may  differ  from 
one  "identical"  equipment  to  another.  Also,  equipments  are  apt  to  age  differ¬ 
ently.  Therefore,  it  is  highly  unlikely  that  the  EMI  performance  curves  of 
"identical"  equipments  will  be  identical.  This  is  illustrated  in  Figure  3 


where  the  EMI  performance  curves  for  many  "identical"  speech  communication 


Figure  3  EMI  performance  curves  for  many  "identical"  speech 
communication  systems  subjected  to  the  same  desired 
and  interfering  signals. 

systems  are  shown.  As  in  Figure  3,  the  selected  performance  criterion  is  I 
articulation  score  which  is  plotted  versus  the  signal-ro-inlerfercnco  utii 
S/I.  In  generating  Figure  3)  it  is  assumed  that  each  equipment  is  exposed 
to  the  same  electromagnetic  environment  (i.e.,  each  equipment  experiences  i 
same  desired  and  interfering  signals).  Observe,  in  Figure  3,  that  at  the 
performance  threshold  the  signal -to-interfcrence  ratio  varies  from  -3'i  dli 


-12  <JB,  depending  upon  the  equipment  under  test.  Consequently,  for  a  given 
desired  sigr.nl,  there  is  a  12  dB  variation  in  the  susceptibility  level,  ].. 
For  example,  when  S  =  -10  dBm,  L  ranges  from  +14  dBm  to  +2  dBm.  In  the 
probabilistic  approach,  the  uncertainty  in  the  susceptibility  level  Is  in¬ 
corporated  by  treating  L  as  a  random  variable. 

In  summary,  the  probabilistic  model  proposed  for  the  F.MC  problem  of 
Figure  2  is  obtained  by  treating 


1)  the  desired  signal,  s(t) , 


as  a  sample  function  from  a  random  pi\ 


2)  Lite  interfering  signal,  i(t),  as  a  sample  function  from  a  random 


cl  ocess  , 

1)  the  equipment  as  a  sample  from  an  ensemble  of  randomly  distribute.! 
"identical"  equipments  whose  EMI  performance  curves  r<ay  be  characterized  in 
terms  of  Lite  random  susceptibility  level,  L.  For  a  specified  desired  siuna 
L  indicates  tin;  minimum  level  of  the  interfering  signal  above  which  KM1  eceu 
The  remainder  of  this  report  is  devoted  to  use  of  the  probabilistic 


model  in  determining  EMC  and  EMI. 
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3.  the  given  sicnai.  and  ixtkkferlnce  i:mi  and  emc  probability  functions 

For  simplicity  in  developing  the  proposed  prob.ihi 1 i st  ic  EMC  model,  dis¬ 
cussion  is  again  centered  around  the  simple  EMC  problem  depicted  in  Figure  2. 
For  ease  of  discussion,  it  is  assumed  that  knowledge  of  only  the  desired  .md 
interfering  signal  average  powers  are  required  to  determine  EMI  pei f ormance . 

Of  course,  depending  upon  the  particular  problem  at  hand,  knowledge  of  other 
signal  pamneters  such  as  peak  amplitude,  energy,  or  peak  power  may  be  needs.!. 

Assume  s(t),  the  desired  signal  at  the  equipment ,  is  known.  l  et  uncer¬ 
tainty  in  the  equipment's  EMI  performance  curve  be  characterised  bv  the  rat..! 
susceptibility  level,  L.  (In  this  report  random  variables  are  denoted  by 
upper  case  letters  while  values  assumed  by  the  random  variables  are  represc.r  v 
by  lower  C3se  letters.)  Let  the  average  powers  of  the  desired  and  interfere', 
signals  be  denoted  by  the  random  variables  S  and  I,  respectively.  Given 
S  “  s,  the  random  variable  L  equals  the  average  power  of  the  interfering  si  a. •' 
required  to  cause  the  equipment  to  operate  at  its  performance  threshold. 

The  random  behavior  of  1,  is  governed  by  its  condiltonal  cumulat ive  dis¬ 
tribution  function  (CDF)  which  is  denoted  by  F  ,  (?|s).  By  definition, 
Fl|s<l|»)  is  the  probability,  given  S  =  s,  that  the  random  variable  I.  is  los, 
than  or  equal  to  the  value  Z.  Symbolically,  this  is  written  as 

FL|g(£|s)  =  P[L  <_  fc | s ] .  (!) 

In  general,  it  is  difficult  to  theoretically  determine  the  conditional  CDF 
of  L.  However,  it  can  be  obtained  empirically  through  laboratory  measurements 

For  this  purpose,  it  is  necessary  to  have  a  large  number,  N ,  of  "identic. i 
equipments  all  exposed  to  the.  same  desired  and  interfering  signals.  Let  N  ^  ( 
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denote  the  number  of  equipments  experiencing  KM  I  when  the  average  power  of  the 
interfering  signal  equals  The  experiment  proceeds  by  1)  holding  the 
desired  signal  average  power  at  the  constant  value  S  =  s,  2)  gradually  in¬ 
creasing  the  average  power  of  the  interfering  signal,  and  3)  recording  the 
value  of  (2).  A  typical  result  of  the  experiment  is  shown  in  Figure  U. 


Figure  4  Experimental  result  in  measurement  of  conditional 
CDF  of  L. 

Note  that  N  ^ j (0)  =  0  because  £  =  0  implies  the  absence  of  an  interfering 
signal.  The  value  of  the  interference  average  power  at  which  the  first 
equipment  from  the  group  of  "identical"  equipments  reaches  its  performance 
threshold  is  given  by  5  =  £^ .  As  the  average  power  of  the  interfering  signal 
is  increased,  more  and  more  er  the  equipments  experience  FM1 .  Finally,  the 
Nth  and  last  equipment,  fiom  the  group  of  "identical"  equipments  reaches  its 
performance  threshold  when  the  interference  average  power  equals  £^ .  Of  course. 
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an  equipment's  susceptibility  level  must  be  less  than  £  when  it  suffers  EMI 
with  the  interference  average  power  equal  to  some  value  £ .  It  follows 
that  is  an  estimate  of  P[L  <  £|sj.  Dividing  the  ordinate  in 

Figure  4  by  N  and  smoothing  the  resulting  curve  yields  an  approximation  for 
the  conditional  CDF  of  L,  as  shown  in  Figure  5. 


Figure  5  Approximation  to  conditional  CDF  of  L. 

The  curve  in  Figure  5  is  conditioned  on  the  presence  of  a  desired  signal 
s(t)  whose  average  power  is  given  by  S  =  s .  In  general,  as  different  sample 
functions  from  the  ensemble  of  s(t)  excite  the  equipment,  the  value  of  S  changes 
and  a  family  of  curves  for  the  CDF  of  L  is  generated.  A  typical  family  is 
shown  in  Figure  6  where,  as  would  be  expected,  the  susceptibility  level  is 
seen  to  increase  with  increasing  values  of  the  desired  signal  average  power. 
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Figure  6  Family  of  conditional  CDF's  of  L. 

The  conditional  CDF  of  L,  FL|S^'S^’  ciin  6<vcn  a  second  interpretation 
Assume  the  average  powers  of  both  the  desired  and  interfering  signals  are  known 
In  particular,  let  S  *  s  and  I  =  i.  EMI  results  when  the  susceptibility  le.vi 
Is  less  than  the  average  power  of  the  interfering  signal.  Note  that 

FL|s(i|s)  =  P[L  <  i|s]  (2) 

is  the  probability  that  the  susceptibility  level  is  less  than  or  equal  to  i, 
given  S  -  s.  Hence,  F^jg(l|s)  equals  the  conditional  probability  of  EMI  undci 
Circumstance  that  S  =  s  and  1  -  i.  Therefore,  FL|g(i!s)  *s  referred  to  as  the 
given  signal  and  interference  EMI  probability  function. 

On  the  other  hand,  EMC  exists  when  the  succcpl ibil itv  level  is  greater 
than  the  average  power  of  the  interfering  signal.  The  given  signal  and  it. ter' 
EMC  probability  function  is  defined  to  be 

Cj|s(i!s)  -  PI  1.  >  i  1  si  -  1  -  P,  |s(i|s). 


/ 


the 


(3) 


By  definition,  C j.  | g ( i | s)  is  the  probability  that  the  susceptibility  level  is 
greater  than  i,  given  S  -  s.  Consequently,  j g C i { s)  equals  the  conditional 

probability  of  EMC  under  the  circumstances  that  S  a  s  and  I  *»  i.  Since 
CTi_(i|s)  is  the  complement  of  a  CDF,  it  is  a  nonincreasing  function  whose 
value  ranges  from  unity  to  zero  with  increasing  i.  A  typical  plot  of  the 
given  signal  and  interference  EMC  probability  function  is  shown  in  Figure  7. 


Figure  7  Typical  plot  of  conditional  EMC  probability  function 
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A .  DERIVATION  OF  EQUIPMENT  EMC  IN  A  RANDOM  ELECTROMAGNETIC  ENVIRONMENT 


The  given  signal  and  interference  EMI  and  EMC  probability  functions  are 
applicable  only  if  the  average  powers  of  the  desired  and  interfering  signals 
are  both  known.  When  the  equipment  operates  in  a  random  electromagnetic  en¬ 
vironment,  S  and/or  I  may  not  be  known.  Useful  results  arc  then  obtained  by 
performing  ensemble  averages  over  the  unknown  random  quantities. 

Assume  the  probability  density  functions  (PDF's)  of  S  and  I  are  denoted 
by  fg(s)  a°d  f^(l),  respectively.  The  condir  -nal  PDF  of  the  susceptibility 
level,  L,  is  obtained  by  differentiating  the  conditional  CDF  of  L  with  respect 


to  H.  In  particular. 


,  dFL|S<*ls> 

fL|S(^s)  ”  d5l 


In  general,  the  conditional  PDF  of  L  varies  from  one  value  of  s  to  another. 

A  typical  family  of  curves  for  fL|s(^ls)  is  shown  in  Figure  8.  It  should  be 


^LIS  ^ 


U 10) 


Figuie  3  Family  of  conditional  PDF’s  of  I.. 


fL|S  UlS|) 
fusWlSg) 

I s3 ' 

I S4) 
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noted  that  L  and  S  are  likely  to  be  statistically  dependent  random  variables 
even  when  S  and  I  are  statistically  independent. 

The  first  case  considered  is  the  situation  in  which  S  =  s  is  known  but 
I  is  unknown.  Averaging  i g ( i | s)  over  all  possible  values  of  I,  there  result 


Cs(s) 


CL|s(i|s)  f I(i)di . 

—00 


(■>) 


Since 


CI|s(i|s)  =  P [L  >  i | s] 


A  fL|sa'S)df’ 


(6) 


it  follows  that  Cg(s)  can  be  expressed  as 


Cg (s)  -  j°°  j  fL|sa|s)  fI(i)dt  di.  (7) 

—03  i 


Given  S  =  s,  Cg(s)  is  the  expected  probability  of  EMC  obtained  by  averaging 
over  all  possible  values  of  the  unknown  interfering  signal  average  power. 
Consequently,  Cg(s)  is  referred  to  as  the  given  signal  EMC  probabil ity 
function.  As  implied  by  the  notation,  Cg(s)  is  a  function  of  the  part  n  il 
value  of  the  known  desired  signal  average  power. 

A  second  case  arises  when  S  is  unknown  but  it  is  know  that  I  -  i . 
Averaging  i g (i | s)  over  all  possible  values  of S,  we  obtain 


Cj.(i)  = 


#0° 

C j | s ( i I s )  fg(s)ds. 

— oo 


(«) 
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Substitution  of  (6)  into  (8)  yields 


cx(i)  = 


i 


J  £l|s^Is^  fs(s)d£  ds 


/CO  ,co 


fL  S(£,s)d£  ds 


i 


(9) 


where  f  (£,s)  is  Che  joint  PDF  of  the  random  variables  L  and  S.  C  (i)  is 

L  |  b  X 

referred  to  as  the  given  interference  EMC  nrohabil  irv  function.  Giver.  I  =  i, 
C^Ci)  is  the  expected  probability  of  EMC  obtained  by  averaging  over  all 
possible  values  of  the  unknown  desired  signal  average  power.  Analogous  to 
Cg(s),  C  (i)  is  a  function  of  the  particular  value  of  the  known  interfering 
signal  average  power . 

Finally,  consider  the  case  in  which  both  S  and  I  are  unknown.  Averaging 
ci,sUl9)  over  all  possible  values  of  S  and  I,  there  re ?u' ' s 

C  *=  I  Cj.  j  s  (i  |  s)  fg  x(s,i)ds  di  (10) 

—CO  _CO 

where  f  (s,i)  is  the  joint  PDF  of  the  random  variables  S  and  I.  Tvpicallv, 

S  and  I  are  statistically  independent.  Then 

fs,i(s,i)  =  fS(s)  fi(i)*  (11) 

Use  of  (6)  and  (11)  in  (10)  yields 
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C 


r°° 


— oo 


fCO 

CI I S  ( 1 1  s>  V*>  di 

— oo 


fL  I  SCP'  I  s)  fS(s)  fi(i)lU  ds  di 


(12) 


r 


r°° 


roo 

fL  S(£,S)  fI(1)  dl  dS  di’ 
i 


C  is  referred  to  as  the  EMC  probability.  It  is  the  expected  probability  of 
EMC  obtained  by  averaging  over  all  possible  values  of  the  unknown  desired 
and  interfering  signal  average  powers.  In  analogy  to  the  concept  of  re¬ 
liability,  C  is  also  referred  to  as  the  system  compatibility. 

By  definition,  the  probability  of  EMI  is  the  complement  of  the  probability 
of  EMC.  It  follows  that 

1  -  Cg(s)  =  the  given  signal  EMI  probability  function, 

1  -  Cj(i)  «  the  given  interference  EMI  probability  function , 

and 

1  -  C  *  the  expected  EMI  probabi 1 itv . 

To  elaborate,  given  S  =  s,  (1  -  Cg(s)]  is  the  expected  probability  of  EMI 
obtained  by  averaging  over  all  possible  values  of  the  unknown  interfering  signal 
average  power.  Similarly,  given  1  =  i,  [1  -  Cx  ( i ) ]  is  the  expected  probability 
of  EMI  obtained  by  averaging  over  all  possible  values  of  the  unknown  desired 
signal  average  power.  Finally,  (1  -  C)  is  the  expected  probability  ot  EMI 
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i 


obtained  by  averaging  over  all  possible  values  of  the  unknown  desired  and 
interfering  signal  average  powers. 
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5.  THE  EMI  RATE  FUNCTION 


The  EMI  rate  f  uiic  t  Ion ,  for  a  specified  desired  signal,  is  denoted  1  y 
X^(i) .  By  definition,  given  S  =  s,  As(i)di  equals  the  probability  an  equ 
ment  will  experience  interference  for  i  <  I  £  i  +  di  provided  compat ib i 1  it 
exists  with  I  i,  To  be  more  specific,  let  and  denote  the  condit  i 
events 

:  equipment  experiences  EMI  with  i  <  I  <_  i  +  di,  given  S 
Aj :  equipment  does  not  experinece  EMI  with  I  i  ,  given  S  =■  s. 
In  terms  of  the  events  and  , 

PfA  A  ) 

As(i)di  =  p(a1|a2)  =  p(A^y-  • 


Since,  given  S  *  S,  the  joint  event  A, A.,  is  equivalent  to  t:<  <  ■  ,;r. : ;  t  '  -■ 
susceptibility  level  being  in  the  interval  (i,  i  +  di),  it  !hv..,  that 

P(A1A  )  =  P[i  s  L  <_  i  +  di  S  -  s)  =  f  ,r(\  Adi. 


In  addition,  given  S  =  s,  the  event  A,  is  equivalent  to 
susceptibility  level  being  greater  than  i.  Theret ore , 

P(a2)  =  Cj  js(i|s)  =  1  -  Fj  s(i  -O  . 
Subsitntion  of  (3A)  and  (15)  into  (11)  yields 


f,  I  (i  •  )di 

A  (i)di  = 

s  1  h  !s 0.*w 
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Dividing  through  by  di,  we  obtain 


Xs(i) 


fL|S(1’S) 

1__  ^L.|S(ilsT 


dFi.  |  s  ( 1 1  s> 


(17) 


Making  the  change  of  variables  i  =  i  in  (17)  and  rearranging  terms,  the 
differential  equation  for  F^  j  ^ ( ? | s)  becomes 


tlFL|S(g!S) 

df. 


A„(«  Fj Js(e|s)  -  Ar«). 


Ob') 


Because  no  interference  results  when  the  interfering  signal  average  power  is 
zero,  the  associated  initial  condition  is 


Fl|s(°U!  -  0. 


(l 


The  resulting  solution  of  the.  differential  equation  for  F  |  fc|s)  is  given  tv 

L )  S 


f« 

r 


Fi(?|s)  =  l-  e  ° 


A  (u)du 
s 


L  S' 


(20  j 


The  complement  of  !'L  |  g ( ^  I p )  is  C'x|s^2''s^'  Hcnce’ 


CI ) S(i ! S)  =  1  ~  FL|S(i|s)  =  6 


A  (u)du 
s 


an 


Observe  that  knowledge  of  the  EMI  rate  function  A  (i)  is  sufficient  to 

s 

completely  determine  the  given  signal  and  interference  EMI  and  EMC  probability 
f  unet  ions . 

The  conditional  PDF  of  L,  fL|s(^ls).  is  also  completely  determined  by 
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the  EMI  rate  functici.  Differentiation  of  (20)  with  respect  to  P.  yields 

dF,  ,c<A|s) 

L|SX'V|“/  d£  Xs(e)  6 


f,  |c(^|s)  = 


X  (u)du 
s 


(22) 


This  result  can  also  be  derived  by  direct  substitution  of  (21)  into  (17). 
Similarly,  substitution  of  (21)  into  (5),  (8)  and  (10)  results  in  the  ex¬ 
pressions 


Cs(s) 


o 


X^ (u) du 


—03 


fjCiJdi, 


(23) 


and 


CjU) 


f 


—  OO 


r1 

X  (u)du 
s 

e  °  fg(s)ds. 


(24) 


— OO  — OO 


X  (u)du 
s 

fg  I(s,i)ds  di. 


(25) 


Assuming  statistical  independence  between  S  ami  I,  it  is  readily  seen  from 

(21 )— (25)  th3t  knowledge  of  X  (i),  f  (s)  ,  and  f  (i)  is  sufficient  to  determine 

S  S  L 

CT |  fi|s)  ,  C„(s),  C  (i),  and  C.  Note  that  X  (i)  characterizes 
the  equipment  while  f^,(s)  and  f^(i)  characterize  the  desired  and  interfering 
signals,  respectively. 

The  concepts  introduced  above  are  now  illustrated  by  two  examples. 

Example  1:  Assume  S  and  I  are  statistically  independent  uniformly  distributed 


l 
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random  variables  with  PDF's 


and 


fs(s) 


0  <  s  <  s 
—  —  o 


otherwise 


fI<i) 


_1_ 

i 

o 

0 


0  <  i  <  i 
—  —  o 


otherwise . 


(26) 


(27) 


Also,  assume  the  EMI.  rate  function  is  given  by 

Xs(i)  =  Ao(l  -f).  U8) 

o 

Observe  that  the  EMI  rate  function  is  a  constant  with  respect  to  i  but  de¬ 
creases  linearly  with  respect  to  s.  It  follows  that 


X  <u)du  =  X  £(1  -  — ).  (29) 

S  O  s 

o  o 


Substitution  of  (29)  into  (20)  and  (22)  yields 


-X  51(1 
o 

e 


and 


-xom  -  f-) 

‘V1 

'  o 


(30) 


(31) 


The  susceptibility  level,  L,  of  the  equipment  is  seen  to  obey  an  exponential 
distribut ion .  Observe  that  Z  varies  from  0  to  even  though  the  range  over 


which  1(0  is  non /civ  tMtiitls 

1  i  v>:n  0  to  i  . 

o 

Kiln  r «  I 

the  given  signal  anvl  inn  *  n  i\*n 

;c t-  KMC  probaKi  1 

itv  lunct 

i  O  -  - 

O  h 

r> 

C,js(i|s) 

-  t.' 

o 

A  sketch  of  C  i  (ijs)  is  shown 
X  |  S 

in  figure  9a. 

In  the  a 

Cr,o  (ilS  )  1 

Co  (si 

Figure  9  Sketches  of  a)  C  ,  (i's),  b )  C„  ( .■>),  ji; .!  c)  (  'i). 

i  S  I 

the  probability  of  F.MC  should  count  unitv.  Note  ih.it  C  i..  .1  - 

I  S  ■ 

pated.  Also,  reflecting  the  fact  that  the  EMI  rat.  futn-t  i  u:  I ; 

s  =  s^,  c  j  |  s  ^ 1 !  So  ^  c1ual:i  llnl,y»  independent  of  i.  given  s  <  ,  : 
of  EMC  decreases  for  increasing  values  of  the  interference  average 
given  signal  EMC  probability  function  is  obtained  by  us.  of  (_>7)  at 
(23) .  Specif ioall v  , 
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cs<s) 


i 

o 

o 


e 


-X  i(l  - 

o 


o 


X  i 
o  o 


_1 

(1 


[1 


-X  i  (1  -  ^-) 

o  o  s 

o 


(33) 


Cg(s)  is  sketched  in  Figure  9b.  Because  the  EMI  rate  function  decreases  linear 1> 
with  respect  to  s,  C^fs)  is  seen  to  be  an  increasing  function  of  s.  When 
s  =  s  ,  C  (s  )  •*  1,  as  was  the  case  with  C  .  (i|s).  The  given  interference 

O  O  o  1  I  J 

EMC  probability  function  is  obtained  by  use  of  (26)  and  (29)  in  (24).  In 
part icular , 


C1(i) 


-X  i(l 

o 

e 


-X  i 

o 


rr  [1  ~ e 

o 


]. 


(34) 


A  sketch  of  C^(i)  appears  in  Figure  9c.  As  would  be  expected,  C  (i)  is  a 

decreasing  function  of  i  with  C^(0)  =  1.  By  coincidence, 

-X  i 
o  o 

cs(0)  ■  w  ■  -4-i —  • 

0  o 


(35) 


Consequently,  the  EMC  probability  averaged  over  all  possible  values  of  I,  given 
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S  »  0,  Is  identical  to  the  EMC  probability  averaged  over  all  possible  values 
of  S,  given  I  *=  i^.  Observe  that  S  =  0  and  I  ®  iQ  represent  worst  case  situations 
with  respect  to  the  desired  and  interfering  signal  average  powers,  respectively. 
Finally,  the  EMC  probability  is  obtained  by  use  of  (26),  (27),  (11),  and  (29) 
in  (25) .  This  yields 


C  « 


-X  i(l  -  ^-) 

o  s  ,  , 

°  —  -p-  ds  di 
s  i 
o  o 


J_  fio  1  -  e 
i  I  X  i 

O  1  O  O 


-X  i 
o 


di  =  — 
o 


CT(i)di 


JL 
6o  J 


-X  i  (1  -  ~) 

s  .  °  °  So  ,  r  s 

o  -  ds  =  —  °  C_  (s)ds . 

o  X  i  (1  -  — )  o  '  o 

O  O  6 


(36) 


Unfortunately,  a  closed  forn  expression  is  not  possible  for  C.  However,  a  value 
for  C  can  be  obtained  using  numerical  integration.  Ir.  fact,  C  is  nothing  mere 
than  l/iQ  times  the  area  under  Cj(i)  or,  equivalently,  l/sQ  times  the  area 
under  Cg(s).  _««««»— 

Example  2:  As  in  example  1,  assume  S  and  I  are  statistically  independent  uni¬ 
formly  distributed  random  variables  with  PDF's  specified  by  (26)  and  (27).  Mow, 
however,  let  the  EMI  rate  function  increase  linearly  with  respect  to  i  and, 
as  in  example  1,  decrease  linearly  with  respect  to  s.  Specifically,  let 

X  (i)  -  X  i(l  -  — ) .  (37) 

s  os 


Observe  that 


X  (u)du  =  -g-  (1. - -)  f  , 

s  2  s 

o  o 


(38) 
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Consequently,  the  CDF  and  PDF  of  L  are  given  by 


fl|s(«|s>  -  i 


T  <*  -  f>'2 

o 


and 


fMs<i|.)  -  v<‘  -  f>° 

1  O 


(39) 


(40) 


These  are  recognized  as  the  CDF  and  PDF,  respectively,  of  the  Rayleigh  distri¬ 
bution.  Following  the  procedure  of  example  1,  it  can  be  shown  that 


C!|S(ilS) 


(41) 


Cs(s) 


m 


/firrx: 


(erf 


(42) 


where 


erf (x) 


<> 


(43) 


CT(i)  =  [1 

1  X  i 


(4  4) 


and 
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Once  again,  n  closed  form  expression  for  C  does  not  exist.  As  in  example  1, 
the  probability  of  KMC  equals  unity  in  the  absence  of  interference  and/or  win  n 
S  =  So  corresponding  to  a  zero  KMi  rate  function.  It  follows  that 

C1|S(0I“>  ■  Cl|s(iK,>  “  suo>  ■  ‘:l(0,  "  '■ 

as  can  be  verified  by  examining  (AT)  -  (44). 

Typically,  the  conditional  CDF  of  the  susceptibility  level,  F 
con  be  determined  only  by  experiment.  f^|^(lls),  the  conditional  I’D!'  oi  I., 
is  then  found  by  numerical  differentiation.  However,  to  facilitate  analyst*  .  ■  * 
is  convenient  to  approximate  the  resulting  PDF  by  a  veil  known  distribution,  la 
many  cases,  analytical  expressions  can  then  be  derived  for  Cjigd's)  and  '  d;. 
The  given  signal  and  interference  EMC  probability  function  and  tb.c  FKI  rate 
function  are  tabulated  in  Appendix  A  for  some  commonly  used  PDF's. 
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6  .  THE  CONFAH  HIE  TTY  MARC.  I  N 

As  mentioned  earlier,  compatibility  exists  provided  the  average  power  of 
the  interfering  signal  is  less  than  the  susceptibility  level  of  the  equipment 
involved.  In  some  cases,  the  probability  of  EMC  is  more  readily  evaluated  by 
introducing  the  random  variable 


Y  =  L  -  I.  (46) 

In  terms  of  Y,  compatibility  exists  provided  Y  _>  0.  Y  also  gives  the  amount 
of  average  power  by  which  the  interfering  signal  can  be  increased  before  EMI 
occurs.  For  this  reason,  Y  is  referred  to  as  the  compatibil ity  margin .  In 
this  section,  expressions  are  derived  for  Cg(s),  C^d).  and  C  in  terms  of  the 
compatibility  margin. 

The  given  signal  and  interference  EMC  probability  function  is  defined  in 
(3)  to  be 

CjjgCils)  =  P[L  >  i's]  =  1  -  FL|s(i|s) 

(47) 

rL|S«!=)«- 


C  |  (i | s)  has  the  following  interpretation.  Let  an  equipment  be  selected  at 

1 1 S 

random  from  an  ensemble  of  "identical"  equipments.  Apply  to  the  equipment  a 
known  desired  signal  with  average  power  S  =  s  and  a  known  interfering  signal 
with  average  power  1  =  1.  equals  the  probability  of  EMC.  Alter¬ 

natively,  assume  the  desired  and  interfering  signals  with  average  powers  S  =  s 
ami  I  =  i,  respectively,  are  applied  to  the  entire  ensemble  of  "identical" 
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equipments,  Cj|<,(i|s)  then  equals  the  fraction  of  equipments  that  are  likely 
to  experience  compatibility. 

The  given  signal  LMC  probability  function  is  defined  in  (5)  to  be 
Cg (s)  «  |°°  C^gUls)  fI(i)di 

— <50 

fCO 

l  J  fLjg(£js)d£]  fj.  (i)di  (48) 

-<»  i 

OO  (CO 

-  fL,s«ls)  fjCDdt  di. 

-oj  i 

The  corresponding  region  of  integrat ion  is  the  shaded  region,  R^,  shown  in 
Figure  10a.  In  note  that  £  ^  i  .  Now  introduce  into  (48)  the  change  of 
variable 

y  -  £  -  i.  (49) 


31 


The  new  region  of  integration,  R2<  in  the  i-y  plane  is  shown  in  Figure  10b. 
The  boundary  i  3  i  In  Figure  10a  corresponds  to  the  boundary  y  =  0  in 
Figure  10b.  Also,  the  points  in  R^  for  which  %  i  correspond  to  those  in 
R2  for  which  y  0.  It  follows  that 


Cs(s) 


f00 

0 


r°° 

f l ( s (y  +  fI(i)di-3dy* 

~  co 


(50) 


Assuming  statistical  independence  between  the  desired  and  interfering  signals, 
it  is  possible  to  write 

f 1 1 S ( 1 1 =  fi(i) -  (51) 

The  expression  for  Cg(s)  may  then  be  written  as 


Vs) 


•on 

0 


[J  fL,s(y  +  ils) 

->Oj 


fl|s(i|s)di]dy. 


(52) 


Now  consider  the  random  variable  Y  =  L  -  I.  If  L  and  I  are  statistically  in¬ 
dependent,  it  can  be  shown  that  the  conditional  density  function  of  Y  is  given 


by  [4] 


fL|s(y  +  i|s)  fi!s(i!s)di. 

—00 


(53) 


Finally,  substitution  of  (53)  into  (52)  results  in 


C 


.(s) 


i*00 

0 


fY|s(y ls>dy* 


(54) 


As  anticipated,  the  probability  of  F.MC  given  S  =  s ,  equals  the  probability  that 
the  compatibility  margin,  Y,  is  greater  than  or  equal  to  zero.  C  (s)  has  the 
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following  Interpretation.  Let  an  equipment  In'  selected  at  random  I i eia  an 
ensemble  of  "Identical"  equipments.  Apply  to  the  equipment  a  known  desired,  si 
with  average  power  S  =  s  and  an  interfering  signal  with  unknown  average  power. 
Cs(s)  equals  the  probability  of  EMC.  Alternatively,  assume  the  desired,  signal 
with  average  power  S  =  s  and  the  interfering,  signal  with  unknown  average  p>»wi  r 
are  applied  to  the  entire  ensemble  of  "identical"  equipments.  C  (;•)  then  equ.i 
on  the  average,  the  fraction  of  equipments  likely  to  experience  compatibility 
when  the  experiment  is  repeated  many  times  with  the  interfering  signal  chosen 
randomly  on  each  repitition  of  the  experiment. 

Example  3:  Let  the  conditional  PDF  of  L  be  Gaussian  with  mean  s  and  vavi.au:.- 
a2.  In  particular,  assume 

•  -Jr-  '  2cJ  •  <’■ 

.2rcL 


Also,  let 


f1(i) 


1 


e 


C' 


Hence,  the  interference  average  power  is  also  assumed  to  be  Gaussian  with  :v.m 
2 

and  variance  Oj.  .  Since  average  powers  arc  nonnegative,  s  and  m  arc 
assumed  to  be  sufficiently  large  such  that  (r>5)  and  (5b)  are  negligibly  sm.i'1 
for  negative  values  of  l  and  i,  respectively.  Direct  use  of  (48)  yields 


Cs(3) 


r 


i  fi ,|s(^S>  di 


(5 


2,  oL  oL 


|  -  3)‘ 

? 

2o, 

e  ^ 


(i  -  m j ) 


2o  T 


d f  dl. 
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I  lu:  above  expression  for  Cg(s)  is  not  readily  evaluated.  However,  assuming 
I  to  be  statistically  independent  of  boLh  S  and  I.,  the  given  signal  FMC 
probability  function  can  be  determined  from  (54)  where  Y  -  1,  -  I  is  the  corns  it  la¬ 
bility  margin .  Since  L  and  I  are  statistically  independent  Gaussian  random 
variables,  it  follows  that  Y  is  a  Gauss  ian  random  variable  with  mean 


nty  =  s  -  m 


(39) 


and  variance 


2  2  2 
°Y  =  °L  +  °I  ' 


(GO) 


Tboiefere,  the  i onditionul  PDF  of  Y  is  given  by 


(y  -  nty)' 


fY|S(yU)  -  - 


—  c 


2o  Y 


>'2v  o 


(■' :  ; 


Substitution  of  (61)  into  (54)  results  in 


(y  -  My)' 


cs(,)  - 


fYlS(y!8)dy  = 


rc 


2o 


/2t T  Oy  J  0 


Y  dy 


((>:*) 


Introducing  the  change  of  variables 


y  -  m 
z  - - * 


( <>  \ ' 


(G2)  becomes 


VB> 


-  f " 


z 

2 


Y  e 


dz  = 


m 

-2. 

rr 


nv 

erf(-:) 
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where  erf (x)  is  defined  in  (43). 


The  given  interference  EMC  probability  1  unction  is  defined  in  (8)  to  be 


C^i)  = 


j  C^Uls)  fs(s)ds 


r® 


J 

—00 


f00 

[  fL|s(l|B)cU]fs(s)ds 

i 


f00 

fL  s(e,S)dl  dS> 
i 


The  corresponding  region  of  integration  is  the  shaded  region,  R^  >  shown  in 
Figure  11.  Interchanging  the  order  of  integration, 


C-^i)  becomes 


Figure  11  Region  of  integration  for  C^(i). 


•CO 

i 


•oo 

r  (£,s)'.ls]<i£. 

_  co 


(65) 


Cj(i)  = 


(66) 


However,  the  marginal  PDF  of  L  is  related  Lo  tiic  joint  PDF  of  I  and  S  according 


to  the  relation 


f,  (£)  =  f  fj  Sa,s)ds.  (67) 


Substitution  of  (67)  into  (66)  yields 


Cj(i) 


f, 
i  ^ 


(-S) 


where  i  is  the  known  average  power  of  the  interfering  signal  .  Now  consider 
the  random  variable 

Y  =  L  -  i  (69) 

where  i  is  recognized  to  be  a  known  constant.  The  conditional  PDF  of  Y  is  re¬ 
lated  to  the  marginal  PDF  of  L  as  follows: 

fYjI(y|i)dy  -  £L(f>)d£.  (70) 

Hence,  the  given  interference  EMC  probability  function  can  be  expressed  as 

6  (i)  =  (  f  i  (yli)dy.  (711 

1  Jo 

Analogous  to  (54),  the  probability  of  EMC,  given  T  -  i,  equals  the  probability 
that  the  compatibility  mat  gin,  Y,  is  greater  than  or  equal  to  zero.  C  (i)  has 
t'’0  following.  Interpretation.  Eel  an  equipment  be  selected  at  random  from  an 
ensemble  of  "Identical"  equipments.  Apply  to  ibc  equipment  a  known  interfering 
signal  with  .average  power  I  -  i  and  a  desired  sigit.il  with  unknown  average  power. 


36 


C  (i)  equa  1  s  the  probahil  itv  of  KMC.  Alternatively,  assume  the  int  i  r 1  .  r  i  a  .• 
signal  with  average  power  I  =  i  and  the  desired  signal  with  unknown  aver*, a  , 
are  applied  to  the  entire  ensemble  of  "identical"  equipments.  C.  (i)  then  < - ‘ ! 1 
on  the  average ,  the  fraction  of  equipments  likely  to  experience  e.  i  :  ; 

when  the  experiment  is  repeated  many  times  with  the  desired  signal  .  ,  , 
randomly  on  each  repitition  of  the  experiment. 

Finally,  assume  1  to  be  statistically  independent  of  both  S  and  .'he 

KMC  probability  is  defined  in  (1?.)  to  be 


j  C'ijs(i!s)  i\.(s)  f  ^  ( i)  ds  di 


-CO  /OO  #co 

=  j  [  fLjs(l|s)di]fi.(s)  f1(i)ds-di 

—CO  —CO  i 


j  fLjS(f.,s)  f I(i)dJ  ds  di. 


—CO  -OD  £ 


Integrating  first  with  respect  to  s,  C  is  given  by 


j  J  [ J  fj  g(J ,s)ds] f x (i)  d •  di . 

—  OO  \  —CO 


Substitution  of  (67)  into  (73)  yields 


fj  (0  fjCDiH.  di. 


Observe  the  simlarity  betv/oen  (74)  and  the  last  equation  in  (48).  'I  lie  onlv 


difference  is  that  the  conditional  l’DF  of  L  appears  in  (48)  whereas  th 


e  mar 1  n  1  I 
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PI' I'  of  L  appears  in  (74)  .  Applying  the  sane  procedure  used  to  convert  (48) 

Into  (54),  it  follows  that 

C  =  |  fy(y)dy  (75) 

0 


where  Y  =  L  -  I.  Once  again,  the  probability  of  EMC  equals  the  probability 
that  the  compatibility  margin,  Y,  is  greater  than  or  equal  to  zero.  Note  that 
C_'s)  in  (54),  CT(i)  in  (71),  and  C  in  (75)  use  the  PDF's  fv|c(>’|s),  f.,|T(y'il, 
and  f  (y)  ,  respectively.  C  has  the  following  interpretation.  Let  an  equipmer. : 
be  selected  at  random  from  an  ensemble  of  "identical"  equipments.  Apply  to 
the  equipment  desired  and  interfering  signals  whose  average  powers  are  unknown. 

C  equals  the  probability  of  EMC.  Alternative! y ,  assume  the  desired  and  inter¬ 
fering  signals  whose  average  powers  are  unknown  are  applied  to  the  entire 
ensemble  ot  "identical"  equipments.  C  then  equals,  on  the  average,  the  fraction 
of  equipments  likely  to  experience  compatibility  when  the  experiment  is  repeat* w 
many  times  with  the  desired  «p d  interfering  signals  chosen  randomly  on  mrii 
repetition  of  the  experiment. 

Example  4:  As  in  example  3,  all  random  variables  are  assumed  to  be  Gaussian. 

In  particular,  “.he  PDF's  of  fT  i_(J.|s),  f.(s),  and  fT(i)  are  approximated  bv 

L 1 5  b  i 

U  -  s)2 


_ ]_ _ 

/2tT 


(76) 


O 


(  7  7 ) 
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(i  -  m  j ) " 


fl(1)  -  * 

'  2 " 


Given  (7o)  ,  (77),  and  (78),  it  Is  desired  to  evaluate  C^(i)  and  C  us i tig  (71) 
and  (73),  respectively. 

The  first  step  is  to  determine  the  marginal  PDF  of  I..  With  reference 


to  (67), 


fLU)  =  fL,S(t’s)ds- 


However , 


fL:s(;ls)  Vs> 


27r  °h  as 


if  (g  -  s)2  "_V_  . 

"  2  a  2  ~c  2'~ 

,  °L  CS 


By  completing  the  square  in  the  exponent  with  respect  to  s,  it  can  be  shown 


2  2  .2  2 

.  ir  °S  +  °l  u  *°S  '  WS  °L  ■■ 

2  2  2  s”  2  2 
.  I  °L  °S  °S  +  °L 

1l,s(£*s)  =  TtToT^  e 


1  (?  -  ms) 

2  2,2 

°S  +  °L 
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Substitution  of  (SI)  into  (79)  results  in 


1 

(f  - 

ms)2 

2 

2 

2 

+  0. 

s 

l 

V 

i  a, 

CJ,. 

L 

,  s 

1 

u  - 

msr 

2 

? 

2 

0„ 

+  0. 

s 

L 

c 

2v 

a.  c>. 

l. 

S 

i 

/nr  /os2  +  oL2 ' 


—CO 


/2ii' 


2 

] 


ds 


(82) 


Therefore,  L  is  a  Gaussian  random  variable  with  mean  and  variance  (c^2  +  C?  Z)  . 

Since  average  powers  are  nonnegative,  m  and  m  are  assumed  to  be  sufficiently 

b  JL 

large  such  that  f^s),  fL  g(£,s),  and  f^(£)  are  negligibly  small  for 

negative  values  of  their  arguments. 


To  evaluate  C^(i)  by  means  of  (71),  it  is  necessary  to  determine  f^|j(y;i) 

where  Y  =  1.  -  i  and  i  is  a  known  constant.  Since  L  is  Gaussian,  the  conditional 

2  2 

TDK  of  Y  is  also  Gaussian  with  mean  (inc  -  i)  and  variance  (o  +  o.  ). 

b  b  Ij 
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Spe  c  if iea] 1 y  , 


fY  1  l<‘y  I  1  >  = 


J  fy  -  (rng  -  i)  }" 

2  2  2 
°S  +°L 


■/2n 


4? 


S  +  CL 


(8  3) 


Applying  rlic  same  procedure  used  to  evaluate  (62),  it  follows  that 


CL(i) 


r 

Jn 


fy | x (y ! i)dy  =  erf ( 


ins  ~  1 

r\  ,  2' 

/oS  +  °L 


(84) 


where  erf(x)  is  defined  in  (43), 

Finally,  to  evaluate  C  by  means  of  (75),  it  is  necessary  to  determine 

f  (y)  where  Y  =  L  -  I.  Let  1  be  statistically  independent  of  L.  Since  L 

and  I  are  both  Gaussian,  the  marginal  PDF  of  Y  is  also  Gaussian  with  mean 

2  2  2 

(m  -  ns  )  and  variance  (oc  +  o  +  oT  ) .  In  particular  , 
b  1  b  L  L 


fy(y) 


^  [y  “  (rng  -  m  )1 

2  2  2  2 
°S  +  °L  +  °I 


/T-  /  2  ,  2“  2 

/2tt  /o„  +o,  +  o. 


(S‘>) 


Applying  the  same  procedure  used  to  evaluate  (62)  and  (84),  it  follows  that 


m  m, 

C  =  j  fy(y)dy  =  erf  (  - ^  ) 

/oS  +  °L  +  °I 


(KM 
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It  is;  of  interest  to  compare  the  expression  for  Cg(s)  1°  example  3  with  these 
obtained  for  C ^ ( i )  .aid  C  in  example  4.  In  each  rase,  the  expressions 
involve  identical  functions  of  the  means  and  variances  of  f y |^ ( v Is >  ’  * v j ]  0' 1  * ) < 
and  fy(v)  ,  respectively. 

It  is  not  always  possible  to  obtain  analytical  expressions  for  C  (s) , 

C  (i)  ,  and  C,  even  when  well  known  PDf-’s  are  specified  for  f<.(s), 

and  f  (i) .  Numerical  integration  techniques  must  then  be  used  to  evaluate  the 
integrals  in  (5),  (8),  and  (12)  or,  equivalent  1 y ,  in  (54),  (71),  and  (731. 

To  illustrate  some  cases  where  analytical  expressions  are  possible. 


cs(s) 


/OO 

cx [s(1 I s)  fi(i)di 

—CO 


fY|s(yis)dy 

0 


(87) 


is  evaluated  in  Appendix  3  for  some  well  known  PDF’s  of  t'T  js(l!s)  and  ^(i). 
Consistent  with  the  property  that  the  susceptibility  level  of  an  equipment  js 
likely  to  increase  as  the  desired  signal  average  power  is  increased,  the 
conditional  PDF's  of  L  are  specified  in  Appendix  B  such  that  larger  suscepti¬ 
bility  levels  are  more  probable  with  larger  values  of  s.  This  results  in 
Cg(s)  being  a  nondecreasing  function  of  s. 

Knowledge  of  fy|g(yls)>  fy|^(y|i)*  and  f^,(y)  can  provide  additional  in¬ 
sight  into  C^(r>) ,  C^(i),  and  C,  respectively.  When  certain  symmetries  in  tin- 
underlying  PDF's  exist,  general  conclusions  can  be  derived  concerning  the 


probability  of  KMC*  For  example,  assume  f^|g(A|s)  and  f^.(i)  are  both  sym¬ 
metrical  functions  about  the  same  mean  m.  Therefore,  for  all  x. 


f,  |o(m  +  xls)  =  fT  |c(n>  “  x!s> 


fj(m  +  x)  =  fj(m  -  x)  . 


From  (53)  conditional  PDF  of  Y  =  L  -  I  is  given  by 


fY|S(y|s)  =  J  fL|S(y  +  lls)  fl|s(i|s)di 

— oo 

•oo 

f  |  (y  +  i|s)  fjdJdi 


where  fj|g(i|s)  has  been  replaced  with  f^(i)  by  assuming  statistical  indeper 
between  S  and  I.  Observe  that 


fY|s(-yls)  “  fL]S(_y  +  ils)  fxO)di. 


Introduction  into  (91)  of  the  change  of  variable 


i  =  m  +  x 


and  utilization  of  the  symmetry  properties  stated  in  (88)  and  (89)  results  ii 


*Y|s(-yis)  =  j  +  x  -  yls>  f1(m  +  *)dx 

— OO 

“  f  f. |c(m  -  X  +  y | s)  fT(m  -  x)dx. 


L 


'.3 


Finally,  with  the  than*''  >->1  variable 


w  -  m  -  x , 


(94) 


(93)  become:. 


f 


Y|S 


C-yU) 


J  fl  :  S^V  +  f  1  ^w^(lw 


f  ,  (v'; 

Y  SV 


)  • 


(93) 


Hence,  fy|s(y<s)  is  even  function  of  y. 


It  foil ows  l ha t 


Cs(»)  - 


eJO 

l 


f y i s ( y ! s ) dv  c  I 


(9b  j 


independent  of  the  particular  PDF’s  for  f^|g(t|s)  and  f  ( L )  as  lonp,  as  the 
conditions  in  (88)  and  (89)  are  satisfied.  Similarly,  it  can  be  shown  that 
CjCi)  =  1/2  provided  f L ( X.)  is  a  symmetric  function  about  the  mean  i.  Also, 
C  =  1/2  when  f^(£)  and  fj(i)  are  both  symmetrical  functions  about  the  same 
mean  ra. 
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7.  Mnmm:  r_r.imM<'0N.r!- 

Thus  far,  discussion  has  been  restricted  to  t  lie  special  case  where  the  i  ‘  ’ 
performance  carve  is  subdivided  into  two  pcrf oruuuce  categories  by  a  sinrle 
perf  ort-iance  threshold,  as  sliouai  in  Figure  1.  Tit  certain  appl  icat  ions  it 
be  desirable  to  subdivide  the  EMI  performance  curve  into  three  or  more  per¬ 
formance  categories.  For  example,  it  may  he  useful  to  categorise  perf errant  e 
as  either  1)  acceptable,  2)  marginal,  or  3)  unacceptable,  as  illustrated  in 
Figure  12.  The  i.MC  per f oiti.-j.icc  threshold  is  that  value  of  the  portorrancv 


Figure  12  EMI  performance  curve  of  a  speech  communication  system  suhdiv iced 
into  three  performance  categories. 


criterion  which  demarcates  tin.'  regions  of  acceptable  and  marginal  performance. 
Similarly,  (he  KM!  performance  thre.sliol d  is  that  value  of  the  performance 
criterion  which  demarcates  the  regions  of  marginal  and  unacceptable  performance. 
As  in  the  single  threshold  case,  EMC  and  EMI  are  said  to  exist  when  the  per¬ 
il,  finance  is  acceptable  and  unacceptable,  respectively.  By  definition,  the  EMC 
scscept ibil  i tv  level ,  1.  ,  is  the  level  of  the  interfering  signal  that  results, 
for  a  specified  desired  signal,  in  operation  at  the  EMC  performance  threshold. 

In  addition,  t  ho  EMI_  sti  scent  ibil  it  v  level ,  L^.,  is  defined  to  be  the  level  of 
the  interfering  signal  that  results,  for  a  specified  desired  signal,  in  operation 
at  the  EM!  performance  threshold. 

In  the  probnbi 1  is!  ic  approach,  the  susceptibility  levels,  and  L^,  are 
treated  as  tandem  variables  whose  behaviors  arc  governed  by  the  conditional  CDF's 


and 


f,(;!s<\:W  "  •'  'c  i  1 

fL  |S(«,|*)  -  HlLj  ^  tj|«l. 


ii.e  corresponding  PDF's  arc 


dFJ  jq(^r|s) 

£r  -(Vis)  Lc|s  c 
rc's  c  ‘  dE - 


and 


1 


(tj |s) 


(97) 


(93) 


(99; 


(100) 
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In  the  ensuing  discussion,  it  is  assumed  that  performance  is  degraded  as  the 
average  power  in  the  interfering  signal  is  increased.  It  follows  that  I.  and 
Lj  obey  the.  inequality 

Lc  <  4.  (101) 

EMC  exists  when  the  susceptibility  level,  4 ,  is  greater  than  the  average 
power  of  the  interfering  signal.  Given  S  =  s  and  I  =  i,  the  conditional 
probability  of  EMC  is  given  by 

ci j s(i I s)  =  P[LC  >  i’  LI  >  i!s] 

=  PtLj.  >  i|(I.c  >  i),s]  P[LC  >  i  |  s  ]  . 

Because  of  the  inequality  in  (101), 

PtLi  >  il(Lc  >  i)*  sl  =  lm 
Therefore,  (102)  simplifies  to 

Sis'1'8)  -rn-c  >  d«]  - 1  -  pLc|S(i|.) 

s) 

Cl|S^ilS^  iS  known  as  the  8lven  signal  and  interference  EMC  probabil ity  function. 

EMI  exists  when  the  susceptibility  level,  I,  ,  is  less  than  the  average 
power  of  the  interfering  signal.  Given  S  =  s  and  I  =  i,  the  conditional 
probability  of  EMI  is  given  by 

4|s(i|s)  =  P[br  <  i,  Lc  <  i  |  s  ] 

=  P(I.C  <  i|(I,j  <  i),s]  I’ll^  <  i|s]  . 


(102) 


(103) 


(104) 


(103) 


hi 


Because  of  the  inequality  in  (101), 


PUC  f  i  |  (Lj  £  i).t 


Therefore,  (105)  simplifies  to 


l^lsUU)  -  H 4  <  1|»1  -  rL  ,5(41.) 


(10?) 


U 1 !  S  ( i  ]  s )  is  known  as  the  g  iven  signal  and  interference  EMI  probahi 1  ilv  fund  ion. 

Marginal  performance  exists  when  the  average  power  of  the  interfering  signal 
is  greater  than  or  equal  to  the  EMC  susceptibility  level,  I.  ,  but  is  less  than 
the  EMI  susceptibility  level,  L ^ .  Given  S  =  s  and  I  =  i,  the  conditional 
probability  of  marginal  performance  is  given  by 

M1|S(i|s)  =  P(LC  <  i,  >  ijs) .  (10b: 

However,  performance  is  either  acceptable,  unacceptable,  or  marginal.  Since 
these  are  disjoint  events. 


ci|s(i,s)  +  ui|s(i!s)  +  Mi|s(iis)  =  K 

Solving  for  an^  utilizing  (104)  and  (!07),  it  follows  that 

Mi|s(1h)  ’  1  -  SpUl")  -  «,|S<‘!») 

■  1  -  <1  -  F,.c|s<i|.»  -  U: 

‘  flc|s(1I»>  - 


(10«) 


If.  -  f.  ,  <  ■  1  =0  J.l  ■  . 

1  "c  O  1 '  ‘ 


-4  8 


•  j  <  i  I  :■ )  is  known  tin1  e  i  von  ..  j  gu.il  .in  1  ijnloij  t- 1 


prob  •!>  i J  i  t  v  i'un._t  iojl . 

The  given  signal  and  interference  KMC ,  I'.Ml  ,  and  m  rgin.il  p«  rs  ..•rr..a:;. 
probability  functions  are  applicable  only  if  the  average  power-  '•  the  d 
and  interfering  signals  are  both  known.  When  S  =  s  is  known  but  !  is  un 
the  appropriate  probabi 1  it ies  are  obtained  by  averaging  over  the  unknown 
average  power  ol  the  interfering  signal.  Let  f1(i)  denote  the  PDF  et  I. 
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given  signal  EMI  prohah  i  1  icy  function 
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fi 


— CO  — oo 


ffLc|S('lF)  -  ri.1|S<-'-s)lrifnd  '!f 


-  given  signal  marginal  per f  onr  nice  grebaC 

[ unc  t i on  . 


i 


When  1  =  i  is  known  but  S  is  unknown,  the  nppropr  i  it  o  pn>b  ah  i  1  it  us  arc 
obtained  by  aver.:;.',  i  up  over  the  unknown  average  power  ot  the  arsirt-J  .sipnai  . 
Let  t'  (s)  denote  tin-  bOF  of  S.  then 


t- 1  (  i )  ils(i's)  f  (S)ds 


i,  ,c(?..|s)  !'  (s)jrf  dr. 

!.  .  S  e  t>  C 
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Uj(i) 


-  j \  on  interference  KMC  piobahi’i  it  y  fun_ct_L 


j  l:2  |S(1  |S)  ES(f0d8 
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-  given  interf  orenco  K.'U  prob.-.b  i  :  : 


V‘>  '  j  »l|S<‘!«>  V*><* 


• .  rv.«<c-c  •>-  -  j  j  >■' 


=  given  inter ferenoc  tnn re, i tu  1  per  foi~’»:itu-c  proL.-ihi 
t  nnc  t.  ion  . 


Ki:i.i!!v,  w!>o:j  both  S  .  i : : » i  1  .in'  unknown,  t  ho  appropriate  prol>;i!>  j  )  i  l  if.'-  an 
obt.t  itu  d  bv  avorv.p.ii!}-,  ovor  both  unknown  nvrrarc  power;; .  Assuming  ;  tat  iatiral 
i iHii’peiK’o'.H'c  between  S  and  1,  tliuR1  results 


r  ^ 
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,cc  f« 
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EMI  probability, 
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Comparison  of  (11 1 )— (1 19)  with  the  corresponding  expressions  previously 
developed  for  the  single  performance  threshold  case  reveals  the  two  sets  of 
expressions  to  be  quite  similar.  It  can  be  shown  that  this  similarity  extend 
to  the  general  case  for  which  (n  -  1)  performance  thresholds  are  used  to  sub¬ 
divide  system  performance  into  n  categories. 
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8.  SAM1M.K  SIZE  CONS 1 DERATIONS  J_N  TI!E  KXITK  1  MENTAL  DETERMINATION  or  F  ( f  ) 


Consider,  once  again,  an  EMI  performance  curve  with  a  single  performance 
threshold  such  as  appears  in  Figure  3.  In  the  probabilistic  approach,  the 
performance  curve  for  an  equipment  is  characterized  in  terms  of  the  random, 
susceptibility  level,  L.  The  random  behavior  of  L  is  governed  by  its  ion. lit: 
CDF,  FL|s(i|s),  which  is  defined  in  (1).  An  experimental  procedure  for  dot «.  r 
mining  F^|g(i|s)  through  laboratory  measurements  was  described  in  Section  3. 
It  was  pointed  out,  for  this  purpose,  that  it  is  necessary  to  have  a  largi 
N,  of  "identical”  equipments  all  exposed  to  the  same  desired  and  inti r: m 
signals.  In  this  section,  a  inethoj  proposed  by  Kolmogorov  and  Smirnov  [h-7 
is  presented  for  determining  a  suitably  large  value  of  N. 

Denote  the  empirically  obtained  conditional  CDF  of  L  by  El' . c  ( 1  j  *) .  Accor 

L  j  b 

to  Kolmogorov  and  Smirnov,  as  long  as  F  i  (Si's)  is  a  continuous  CDF,  then 


-  6  <  FL|s(t|s)  <  Ej|s(t|a)  ♦  4]  -  r 


provided 


where  6  <  0.2  and  a  is  determined  from  Table  1  [8] 


TABLE  1 

PARAMETERS  USED  IN  KOEMOCOROV-SM I RNOV  TEST 


0.85 

1.1  A 


0.00 

1.2? 


0.95 


0.99 


1.30 


1  .63 


).  ,i:  a  ici:  (.1.0:  mean  th.it ,  with  j>  rob.ih  !  1  i  I  v  ,  the  itu  411.1 1  i  l  y  ,  I  ^  |  |  fi)  - 

.  ( V.  |  :•)  i  is  su  t  i  • . !  ied  fir  .ill  In  the  iiUitv.il  (-<v,  -*-■)•  In  other  wore 


1  . : 


vli.i t ever  tin-  Cl '1'  r. 


, ,  .,0  'or  is  it  i-.  cent  ii.noub ,  the  random  dmi-iin  b<  - 


tween  the  boundar  ies  l-.j' ,  •  I;.)  -  5  n:u!  F.^  1  ^  (  7  J  s )  4  d  contains  Fj  |<.(ijs)  comple 

wiiii  prohahi  !  i  ly  .  Not'1  th.u  the  noma  in  quoted  depends  on  the  empirical 

distribution  !..  j  ^.  (  >.  |  s )  obtained  using  t  lie  particular  sample  of  N  equipments. 

Should  the  experiment  he  repeated  using  a  different;  sample  of  N  equipments, 

E‘  ,  (il.u)  i  ■  likely  to  change  causing  the  domain  of  width  2  6  to  be  centered 

1,  |  S 

about  a  different  1  unction. 

N  1  N  , 

In  the  statist  ic..  i  literature,  the  interval  (Ej  |g(l|s)  -  6,  1  c  (  ,  s)  f 

is  cal  leu.  the.  c-'hf  id,  n.  e  iiiferv.il  for  estimating  F  j  (f  js)  .  The  boundaries 
E(  jj,(t|s)  -  ‘  and  Ej  jg(«  j  n)  +  f.  are  called  confidence  limits  for  Fj  ^  j  s )  . 

The  probability y  is  c.'i’lod  the  conf  idence  c  oe  f  f  i  c  i  en  t .  Multiplying  ,  bv 


1 0U'  ,  1  (HI  i  1  j.s  called  t  he  con  fi  donee  level .  Thus,  (120)  specifies  a  coniidi.ie 
domain  for  ,n;  unknown  cumulative  distribution  function. 

1  ■.  ciple _ a :  Let  t  =  0 . 1  5  and  y  =  0.90.  What  should  be  the  sample  size,  \' ,  of 

tlic  experiment? 

From  Table  1 ,  a  =  1.22.  Substituting  into  (J21),  there  results 


1  22  ? 

N  -  W  “  66-15-  ( 

:!  in  e  ,  <i  7  equipments  are  needed  if,  with  prohahi  1  ity  0.90,  F^  ,  0  (  •  ,  h  is  t.  f- 
co.  ip  1  e  t  e  I  y  contained  within  the  confidence  interval  of  width.  0.3  centered  on  i 
c.np  i  1  *  cal  CL1'. 

F.  n.qile  L;  Let  ‘  =  0.1  and  y  -  0.975.  What.  should  be  the  sample  size,  N,  of 


is-.pei  fluent  ? 


Note  that  Y  =  0.975  does  nut  appear  in  Table  1. 


i  i  ■  r i  •  o(jii<  ■■ ,  i  1 


it  is 


necessary  to  interpolate  between  the  tabulated  values.  Using.  linear 


0.975  -  0.95 
a  -  1.36 


0.99  -  0.95 
1.63  -  1 .36  ' 


Solution  of  (123)  for  a  yields 


a 


1.36 


(0.025) (0.27) 
(0.04) 


1.53- 


Using  (121),  the  inequality  on  N  becomes 


1  53  2 

N  I  =  23A- 

Therefore,  234  equipments  are  needed  if,  with  probability  C.975,  F  , 

L  |  S 

to  be  completely  contained  within  the  confidence  interval  of  width  0. 
on  the  empirical  CDF. 


i 1 1  :  p<j  la:  i<  r.  , 

(123) 

(124) 

(125) 

( ;  !  s)  is 
centered 


9 .  APPROXIMATION  OP  FL  j  ,,  (  2  |  s)  li V  A  WF1.1.  KNOWN  PI  STR 1  Cl'T  J  ON 


When  j g ^ ^  1  ^ ^  ■*-s  determined  in  the  laboratory,  it  is  convenient,  for 

analytical  reasons,  to  approximate  the  experimental ly  obtained  CDF  by  a  weiJ 
known  distribution.  In  general,  several  different  distribu'.  ions  may  provide 
acceptable  approximations  to  the  experimental  data.  This  section  presents  some 
statistical  techniques  for  rejecting  candidate  dist.ribut ions  which  arc  not 
supported  by  the  data. 

The  procedure  consists  of  two  steps.  First,  two  measures  which  contain 
information  relative  to  the  shape  of  a  distribution,  the  coefficients  of 
skewness  and  kurtesis,  are  used  to  make  preliminary  selections  of  candidate 
distributions.  The  candidates  ate  then  tested  to  see  whether  any  shot! Id  be 
rejected  on  the  basis  of  a  significant  statistical  deviation  between  the 
experimental  data  and  the  distribution  being  tested.  The  approximation  to 
Fl,s<*|s)  is  chosen  frv.m  those  distributions  which  passed  the  test.  However, 
it  is  not  possible  to  st.ate  whether  any  one  distribution  which  passed  the  test 
is  better  than  any  other. 

The  coefficients  of  skewness  and  kurtosis  are  related  to  the  kth  centra! 
moment  of  a  random  variable  where  k  -  •  -.  hot  the  mean  of  the  random 

variabLe  1.  be  denoted  by  m^.  The  k— -  cent  ra  1  moment  of  I.  is  defined  to  he 

Hk  »  E[  (L  -  ntj  )k]  (  i  . 

where  K[  ]  denotes  the  statistical  operation  of  exported  value.  I’.v  definition, 

u3 

a.  r-  — r -  coef  I  icieni.  of  skewness  (  | 

H 


56 


and 


u 

4 

'<*  2 
11 2 


-  coefficient  of  knrlosi: 


The  coefficient  of  skewness  is  a  measure  of  the  asymmetry  of  the  PDF.  The 
coefficient  of  kurlosis  is  a  measure  of  the  "peakednoss"  or  "f  1  n !  in-.ss"  of 
a  PDF  in  the  central  part  of  the  distribution. 

When  experimentally  determining  F^|g(i|s),  it  is  assumed  that  each  of 
"identical"  equipments  is  exposed  to  the  same  desired  and  interfering  signi 
As  explained  in  Section  3,  the  experiment  proceeds  by  i)  holding  the  desir 
signal  average  power  at  the  constant  value  S  =  s,  2)  gradually  inen-ns  is 
average  power  of  the  interfering  signal,  and  3)  recording  tiie  value  of  1; 
where  N  j(t)  denotes  the  number  of  equipments  experiencing  EMI  with  the  sv 
power  of  the  interfering  signal  equal  to  1.  A  typical  result  of  tin  exptri 
is  shown  in  Figure  4.  Let  i,  denote  the  value  of  the  interference  a vein go 
power  at  which  the  i— - -  equipi.it. .t  reaches  its  performance  threshold.  An  uub 
estimator  of  the  mean  m^  is  given  by 


1  N 

“l  =  N  Sl 


i=l 


Similarly,  unbiased  estimators  of  u_,  y  ,  and  y  are 

2  3  4 


1  2 
]12  =  N  -  1  ~  mL)  ’ 


i=l 


_  N 

f'J  r»  —  ^ 

U3  =  (N  -  1)(N  -  2)  i;]<''  i  ~  ’V  ’ 


_ N  r  ..  6N  -9 

•;  '  L  mi  2 

(N  -  1)  (N“  -  3N  +  3)  i~l  N  -  3N  +  3 


"  • 


0  J3) 


Using  (130)-(L32)  in  (127)  and  (128),  estimates  o£  the  coefficients  of  skewness 
and  kurtosis  become 


(133) 


a4  r  ,2 

(»‘2) 


Preliminary  selections  of  some  well  known  distributions  for  approximating 

Fj  |g(fc|s)  are  by  comparing  the  estimates,  e.^  and  u^,  with  the  known  values, 

a.,  and  a  ,  of  the  various  distributions  under  consideration.  Those  disiribut ions 
3  4 

for  which  a  reasonable  match  exists  become  candidates  for  further  consideration. 
In  general,  the  PDF's  of  the  well  known  distributions  contain  two  or  more  para¬ 
meters  for  which  numerical  values  must  be  determined  in  order  for  the  distr  iber  i  i 
to  be  completely  specified.  These  parameters  are  assumed  to  be  chose  i  such  that 

the  mean  and  variance  of  the  distribution  equal,  respectively,  t  lie  sample  mean, 

_  _  2  — 

n>l  ,  and  the  sample  variance,  (a  )  =  .  Therefore,  with  respect  to  the  mean 


and  variance,  there  is  no  loss  in  generality  by  considering  the  standardly 


cu  random 


var iabl e 


(137) 


whore  the  vari.i'ii  c  of  I.  i:'  denoted  hv 


.  Tin  s  t  anda  i  d  i /rd  r.ir.. .00  Vile..'1' 


X,  has  zero  mean  and  unit  variance.  It  can  he  shown  that  Lin-  cod  i  u  .eat'. 


skewness  and  kurtosis  are  invariant  under  standard  i  zat  ion  of  the  rami,):,  varia:  le 


In  other  words,  the  distributions  for  !  and  X,  respectively  ,  have  tin  sa::.e 


numerical  values  of  ct,  and  ct..  For  the  purpose  of  cor.tpar in.:  the  sh--.i  es  of  the 

3  A 


various  distributions,  it  is  convenient  to  work  with  the  st andardized  render’ 


variables.  Analytical  expressions  for  PDF's  whose  random  variables  have  been 


standardized  to  zero  mean  and  unit  variance  are  tabulated  in  Appendix  C.  .'lent. 


with  their  coefficients  of  skewness  and  kurtosis. 


Having  selected  a  group  of  possible  distribut  ions  for  app roxima.t  in: . 


F  I  (0.1s)  on  the  basis  of  comparing  ft,  and  a  with  ft.  and  ,  respective1,  v , 
Lb  j  4  J  q 


it  is  then  desirable  to  perform  a  "goodness  of  fit"  test  on  each  candidate 


distribution.  This  enables  a  distribution  to  bo  rejected  when  j  significant 


statistical  deviation  exists  between  the  experimental  data  and  the  dist r i l  et  ion 


being  tested.  Even  though  a  distribution  is  not  rejected,  It  cannot  be  accepted 


with  confidence.  In  addition,  if  several  distributions  are  not  rejected,  the 


goodness  of  fit  test  cannot  be  used  to  accept  one  distribution  over  another. 


The  best  that  can  be  said  regarding  distributions  which  pass  the  test  is  that 


they  were  not  rejected.  The  goodness  of  fit  test  prevents  acceptance  of  a  d  i  - t t 


bution  which  is  likely  to  make  a  poor  approximation  to  F^  ,  (ijs). 


N,  the  number  of  "identical"  equipments  employed  jn  the  experiment ,  deU-rmi- 


which  goodness  of  fit  test  to  use.  For  small  values  of  N  (i.c., 


30).  the 


Kolmogorov-f mi rnov  test  is  suggested.  On  the  other  hand,  the  "  test  is 


preferred  when  N  is  large  ([.<■.,  N  :•  30) 


Case  1  -  Oral1  N  0»  30):  Once  again,  it.  is  convenient  to  work  with  st  uiclat.. 


1 


b<) 


distributions  having  zero  mean  and  unit  variance.  With  respect  to  the  data,  the 

_  2 

sample  mean,  m  ,  is  given  by  (129)  while  the  sample  variance,  (o^  )  =  *-s 

given  by  (130).  The  experimental  data  is  standardized  according  to  the  relation 

~  ®L 

x  =  - - -  ;  i  =  1,  2,  ....  N.  (136) 

°L 

Denote  the  CDF  of  the  standardized  distribution  to  be  tested  by  F^(x) .  Analytical 
expressions  for  various  standardized  CDF's  are  tabulated  in  Appendix  D. 

The  Kolmogorov-Sini  rnov  test  is  presented  here  without  proof.  It  proceeds; 
as  follows: 

1)  Standardize  the  experimental  data  by  utilizing  (136). 

2)  Obtain  annnalytieal  expression  Fv(x)  for  the  distribution  under  test, 
(Appendix  D  may  be  of  some  value  in  this  regard.) 

3)  For  each  experimental  data  point,  evaluate  FY(x  ). 

A  X 

A)  Construct  an  experimental  standardized  CDF  by  usi  the  relation 


P(XE  <  .p  - 


wy 

N 


(137) 


where  )  denotes  the  number  of  equipments  experiencing  EMT  with  the  average 

power  of  the  interfering  signal  equal  to  a  value  just  slightly  larger  than  f.  .  . 

i 

The  random  variable  denotes  the  experimental  data.  The  experimental  CDF  has 
the  shape  of  a  staircase  function  such  as  shown  in  Figure  A. 

5)  Find  the  maximum  distance  given  by 


D 

max 


max 

i 


W 


•wy 


0  38) 
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6)  Assign  a  numerical  value  to  l  lie  rauilitiuii.il  pi  eb.tb  i  1  i  I  v  whir  h 
is  defined  to  be 

a  =  P  (rejecting;  d istr ibut 1 on | d istr Lbu t ion  is  a  suitable  appro:-:  i  ail  i ..  a  i  .  < 

i 

This  probability  is  known  ns  the  si gnii icanco  1  eve!  of  the  test. 

7)  Refer  to  the  Table  in  Appendix  E.  Find  the  tabulated  value  c-*rr< 

to  the  applicable  sample  size,  N,  and  significance  level,  a.  Denote  the  t  ■. :  - 

value  by  D^. 

N 

8)  Reject  the  distribution  if  D  >  pf . 

max  a 

Example  7:  Assume  10  "identical"  equipments  are  used  in  the  ox  per  it.. cut  foi 
determining  F^|g(i|s).  The  results  of  the  experiment  are  suinmat  i  ::ed  in  i  .0  I.  A. 
The  sample  mean  is  given  by 


"L 


1  10 

J_  V  t 

10  i=l  1 


101 1.1  mw. 


(1  ■”) 


The  sample  standard  deviation  is  given  by 


oL  =  [|  l  (^s.  ~  mL) 21 3 72  =  68 -48  r'lW- 

i=l 

Therefore,  the  standardized  data  are  obtained  using 

li  -  1011. 1 
Xi  =  68.48  ' 

,  3  4 

is  tabulated  in  Table  A  along  with  and  x^  which  are  needed  in  cor; rt 

the  coefficients  of  skewness  and  kurtosis.  Since  x^  has  zero  mean  and  unit 


variance,  it 
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I  V81  !  A 


i  i 

F  ( x  .  ) 

X  l 

p(xk‘V  !fx 

(x.  }-?(  >:  •  X 

1  ! 

.Ai/l 

.3287 

.4 

.0/13 

...  .  .  i.l.'  , 

.6466 

.6 

.0466 

.  >  . 

i! .  i  ;  ■  .  »■>?'.> 

.7774 

.8 

.02  76 

*  >  0  c 

b  .1  ... 

.0723 

.1 

.02  7  7 

,  '(A? 

•ii ,  -  ,;h‘S ;  . O'.ii 

.4  i  23 

.5 

.0873 

0  i  c'JA 

l  .  •  '<  2.4-4  4  2.393 

.8979 

1  .0 

.  1  02 1 

7  9  7  1 

.  |  /]  1  .0936 

.290 

.3 

.010 

8  L0V.I 

.  ;4r-  .  2 187 

.7488 

.7 

.0488 

9  y;j9 

-1.491  A  .99  2 

.0744 

.2 

.1236 

Id  j  OF’  3 

1.0  70  1.138  1.2  Lb 

.8498 

.9 

.0302 

follow,,  that 

O'-' 

ii 

rn 

Is 

\°  ^ 
l  xi 
i=l 

=  -0.406 

( 1  V> ) 

and 

10 

o  =  .ilooi_  y 

4  (9)  (73) 

4 
x . 

1 

,(.5LI  =  1.532, 

(73) 

(  1  M  , 

In  Appendix 

the  Wcilnill  distribution  is  the 

only  distribut ion  for 

which 

nepativo  wilt 

ii’S  of  .'i  appear.  Likely 

candidates  are  those  special  ca 

si's  lor 

which  (•  =  4, 

5,  6,  7,  and  8. 

let  us  Lest  Llie  Vic  i  bull  distribution  for  which  8  =  4.  From  Appendix  ,’i 
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it  is  :CC1)  tll.U 


(x  t  T  .  )  * 

■  239.U 

F^(x)  --  |  1  -  e  )  u(x  +  3.564). 


(145) 


Val  ucs  of  Fv(x.)  ,  !’  (X  x.),  and  |F„(x.)  -  P  (X„  <  x.)  1  are  also  tabulated  in 
X  i  h  —  l  1  X  i  L  —  l 

Table  A.  Observe  that 

D  =  0.1  256.  (U6) 

max 

Assume  the  sigr.il  ic.ince  level  is  set  at  a  =  .05.  Referring  to  Appendix  E, 
it  is  found  that 

=0.409.  (147) 

Since  D  <  D’^  ,  the  distribution  is  not  rejected, 
max  10 

^ ■— ■  ■  ■  * 

> 

Case  2 -barge  N  (N  >  30):  For  large  vale.es  of  N,  the  test  is  preferable.  It 

2 

is  presented  here  without  proof.  The  \  test  proceeds  as  follows: 

1)  Standardize  the  experimental  data  by  utilizing  (1^1 

2)  For  the  distribution  under  test,  divide  the  range  of  X  into  10  cells 
such  that  it  is  equally  probable  X  will  fall  into  each  cell.  In  the  statistical 
literature,  eacli  cell  is  referred  to  as  a  dee  i !  o. .  The  decile:  hem:  „  ir  it  s  for 
some  standardized  CDF’s  are  tabulated  in  Appendix  F. 

3)  Count  the  number  of  standardized  data,  x^,  contained  within  each  cell. 
Denote  the  number  of  data  in  the  k—  cell  by  ;  k  =  1,  2,  ...,  10. 

4)  Compute  the  quantity  1’^  defined  by 


•  •  ) 


10. 


(14b) 
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5)  Evaluate  the  statistic 


10 

l 


k=l 


(Pk  ~  0-D2 

0.1 


(149) 


6'  Let  M  denote  the  number  of  parameters  in  F  (x)  ,  the  CDF  under  test,  whose 

A 

numerical  values  must  be  determined  in  order  for  the  distribution  to  be  completely 
specified.  The  degrees  of  freedom  are  given  by?=10-l-M=9-M. 

7)  Assign  a  numerical  value  to  the  significance  level,  ~t ,  defined  by  (139'. 

8)  Refer  to  the  Table  in  Appendix  C.  Find  the  tabulated  value  corres¬ 
ponding  to  the  applicable  degrees  of  freedom,  P,  and  significance  level,  ct . 

2  a 

Denote  the  tabulated  value  by  (\  )g. 

2  2  a 

9)  Reject  the  distribution  if  v  >  (x")p* 

The  above  discussion  is  applicable  when  N  >  100.  Should  30  <  N  <_  100, 
step  2  should  be  modified  such  that  the  range  of  X  is  subdivided  into  5  equally 
probable  cells.  The  cell  boundaries  can  be  obtained  from  Appendix  F  by  select  in.: 
the  5  intervals  (xQ,  x2>,  (x2>  x^) ,  (x^,  xg)  ,  (xg,  xg) ,  and  (xg,  x1Q) .  The 
statistic  equivalent  to  that  in  step  5  is  now  given  by 


9 

*  N 


5 


l 

k- 1 


( 1  50) 


Finally,  the  degrees  of  freedom  specified  in  step  6  are  now  given  by 
P*=5-l-M  =  4-M.  All  other  steps  remain  the  same. 

Independent  of  whether  10  or  5  cells  arc  employed,  best  results  are 
obtained  when  the  number  of  data  points  in  each  cell  exceeds  three. 


w 


4 


<•  }  :  Au.-unc  20  "  i  ‘.cut  i'\il  ”  equipments  ai'i-  usi  d  in  the  exper  iment  for 

determining,  1'  1 1. I  •  s)  .  lor  th  L:;  value  of  N,  the  Ko  lmugorov-Snirnov  t  e:;t  shut)  1  • 

r> 

be  used .  However  >  the  :•  Lest  will  be  applied.  The  small  value  ol  k  ha..  be<  a 

specified  in  order  to  make  the  example  more  tractable.  The  results  of  the 

experiment  are  summarized  in  Table  B.  The  sample  mean  is  riven  by 


20 


\  20 

i=l 


1007.3j  ir.w . 


The  sample  standard  deviation  is  given  by 


20 


°L  ‘  ll9  ^  ('i  "  V‘ 


=  49.51  in 


Therefore,  the  standardized  data  are  obtained  by  usint 


H  -  1007.35 
Xi  =  49.51 


(: 


n 


( 1  5 


3  4 

x.  is  tabulated  in  Table  is  along  with  x.  and  x  which  ere  needed  in  eonput  i~. 

i  li  r 


the  coefficients  of  skevmess  and  kurtosis.  Since  x^  has  zero  mean  ird  unit 


variance,  it 


TAW  >.  K 


1 

£  (row) 

*1 

y 

‘i 

1 

966 

-.835 

-.582 

.486 

2 

976 

-.633 

-.254 

.161 

3 

1023 

.316 

.032 

.010 

4 

1038 

.619 

.237 

.147 

5 

1064 

1.144 

1.497 

1.713 

6 

1036 

.579 

.194 

.112 

7 

1035 

.558 

.174 

.097 

8 

932 

-1.522 

-3.526 

5.366 

9 

910 

-1 .966 

-7.599 

14.939 

10 

1.033 

.558 

.174 

.09  7 

11 

978 

-.593 

-.209 

.124 

12 

950 

-1.158 

-1.553 

1.798 

13 

1013 

.114 

.002 

.000 

14 

964 

-.876 

-.672 

.  589 

15 

1091 

1.690 

4.827 

8.157 

16 

1058 

1 .023 

1.071 

1 .095 

17 

104  5 

.760 

.4  39 

.334 

18 

1043 

.720 

.373 

.269 

19 

1036 

.579 

.194 

.112 

20 


954 


-1 .078 


-1 .253 


1.350 
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follows  that 


a3 


(20) 

(19) (18) 


20 

l 


i=l 


-.377 


(154) 


and 


(400)  2r°  4  _  (111) 

(19) (343)  i^1  *i  343 


1.945. 


(155) 


With  reference  to  Appendix  C,  it  would  be  logical  to  select  the  Weibull  dis¬ 
tribution.  Howeveri  let  us  test  the  triangular  distribution.  Because  of  the 
small  value  of  N,  let  the  range  of  X  be  subdivided  into  5  cells.  Referring  to 
Appendix  F,  the  5  cells  are  given  by  (-2 .45,- . 900) ,(-. 900,- .25 '),(-. 259, . 259)  , 

(.259,. 900),  and  (.900,2.45).  The  standardized  data  in  each  cell,  for  each  cell, 
2  2 

P^and  (P^  -  0.2)  are  tabulated  in  Table  C.  The  X  statistic  becomes. 


TABLE  C 


Cell  Number 

1 

2 

3 

4 

2 

Cell  intervals 

(-2. 45, -.900) 

(-.900, -.259) 

(-.259, .259) 

(.259, .900) 

(.900,  2 

standardized 
data  in  each 

8,9,12,20 

1,2,11,14 

13 

3,4,6, 7 
10,17,18,19 

5,15,16 

-ell 

\ 

4 

4 

1 

8 

3 

Pk 

0.2 

0.2 

.05 

0.4 

0.15 

(Pk  -  0.2)2 

0 

0 

.0225 

.04 

.0025 

according  to  (150), 


0.2 


6.53 


(156) 


A 
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Since  the  triangular  distribution  contains  two  parameters,  the  degrees  of 
freedom  are  given  by 

P  «=  4  -  M  «=  2 .  (157) 

Assume  the  significance  level  is  set  at  a  '  0.05.  Referring  to  Appendix  C, 
it  is  found  that 

2  0.05 

(X  )2  =  5.99.  (158) 

2  2  0.05 

Since  X  >  (x  ) 2  •  the  distribution  is  rejected. 
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10.  CONCESSION  AM)  Sl'CCESTiON’S  KOk  EtMTHli  MWK 

An  inherent  feature  of  many  IMC  problems  is  the  rai  lomne  s  assoc  iut  ed 
with  the  desired  signals,  interfering  signals,  and  equipments  involved.  In 
this  work  a  probabilistic  approach  is  proposed  which  introduces  the  follow; 
concepts : 

1)  Fj  |  s ( i j  s) ,  the  given  signal  and  interference  EMI  probability 
f unc  t ion , 

2)  C  iQ(i|s),  the  given  signal  and  interference  EMC  probability 

■  I  a 

func  t ion , 

3)  C  (s)  ,  the  given  signal  EMC  probability  function, 

b 

4)  Cj(i),  the  given  interference  EMC  probability  function, 

5)  C,  the  EMC  probability  (also  referred  to  as  the  compatibility), 

6)  X  (i),  the  EMI  rate  function 

7)  Y,  the  compatibility  margin. 

Several  examples  illustrating  the  interrelationships  between  these  concepts 
r re  presented. 

For  problems  where  a  deterministic  approach  is  inappropriate,  the 
probabilistic  approach  can  provide  an  improved  EMC  model.  In  addition,  hv 
utilizing  statistical  parameters,  statistical  maeroimdels  can  ho  developed 
which  simplify  the  characterization  of  complex  signals  and  ct|u  i  pr-ent  .  ih; 
approach  appears  to  be  especially  well  suited  for  the  EMI  mode!  in/,  o  :■  I  r< 
circuits.  Because  of  their  tremendous  complexity,  it  is  noithet  <’,»•,  iw.'n!, 
possible  to  carry  out  an  exact  circuit,  modeling  and  analysis  in  order  to 
determine  EMC. 

Central  to  the  probabilistic  approach  is  the  experimental  dot  i-rr.i  n.i  t  i . 
K|  .(CIS),  A  method  for  determining  a  suitably  large  value  oi  y,  the  imi-Im 


of  equipments  to  be  used  in  the  experiment,  if;  presented.  For  nnalyii.  .< 

reasons,  it  is  desirable  to  approximate  the  experimentally  deletr-ined 

F  | _(" j  s)  bv  a  wol 1  known  distribution.  Some  statistical  techuM  to  ,  :  ■ 

F  |  «•> 

reject  ini;  candidate  distributions  not  supported  by  tiie  experiment  il  data 
are  discussed. 

it  is  recommended  titat  randomness  associated  with  microcircuit  ruse 
tibility  be  invest  ip, a  ted.  To  what  extent  is  Lite  susceptibility  level,  !. 
random  ?  If  the  susceptibility  level  is  random,  can  its  PDF  be  determii; 
Assuming  the  PDF  is  known,  can  the  probability  of  F_MC  and/or  MM  I  be  succ 
fully  determined.’  To  de'  el  op  answers  to  these  questions,  it  is  sugpes:  e 
the  random  susceptibility  of  the  7400  Til.  NANI)  gate,  a  simple  integrated 
circuit,  be  studied  experimentally.  If  the  study  is  concluded  succ essfu 
more  complicated  microcircuits  should  then  be  investigated. 
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Special  Case  of  the  Beta  Distribution 


9.1  a  =  b  =  0  results  in  the  uniform  distribution 


fL|S(£ls)  “ 


2k 

0 


c  <_  £  <  d 
otherwise 
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APPENDIX  B 


EXAMPLES  OF  ANALYTICAL  EXPRESSIONS  FOR  C  (s) 


Given  fg|s(Ms)  <!tut  f  ( i )  ,  the  given  signal  EMC  probability  function  may  E.- 
evaluated  using 


Cs(s) 


cI,s(i|s>  ^(Ddi 


■  r 


where  Cjig(i|s)  and  f  y  |  S  ^ y  ! S  >  are  ^e£<:Jrnn  ned  from  (3)  and  (33),  respective!; 


1.  f1|s(£|s)  Uniform  and  fj(D  Uniform 


r 


fL|Sai»>  = 


i 

2a 


s  -  a  <  L ,  <  s  +  a 


0  ,  otherwise. 


fj(i>  = 


2b 


in  --  b  <  j  <  m  +  h 

otherwise . 


(ijs)  Weibull  and  f^(i)  I’niform 


(tis)  =  |  (t  -  s)K  '  1 


(*  -  s)K 

a 

u(lt  -  s) 


m-b<^i<ra  +  b 
otherwise 


1/K  ,  ,.n  ^  ,  nX  +  1  .  rJC+1 

a _  r  (~i)  r , g  +  b  -  s.  .m  -  b  -  ss  , 

2b  L.  n!(nK  +  1)  U  1/K  1  1/K  ;  1  ’  b 

n»0  a  a 


,  ,  i/K  °° 

s  +  b  -  m  ^  a _  y 

2b  2b’  nKnK  +  1)  v  1/K~7 

ri"U  a 


(-D 


m  +  fc-svnK+l 

' - T777 — )  »  HI  -  b  •:  3  <■  D  +  b 


1 


tn  +  b  <  s 
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APPENDIX  C 


M.ALYTICAL  EXPRESSIONS  FOR  STANDARDIZED  PDF's  AND  THEIR  COEFFICIENTS  OF  SKEWNESS 


AND  KURTOSIS 


Uniform 


ANALYTICAL  EXPRESSION 


,  -/3  <  x  < 


O  ,  <  I : 


Gaussian 


f  +  TS*  1  x  <  o 


Triangular  fx(x)  «j-|  +  0  £  x  <  /E 


0  ,  /S  < 


f  (x)  -  - ± -  e 

C/2?(x  -  b) 


In(x-b) -Inal 


u(x-b) 


2  2 

Log  Normal  a  -  [eC  (eC  -  l)]-*^ 


c2  ,wc2 


(e  -  1)  (e  +  2)  ( e  +  2e 


4C2,  „  3C2 


b  -  -(eC  -  l)'l/2 


2  r 

+3e  -3) 


Special  Cases 


Log  Normal 


C  -  1/4 
C  -  1/2 
C  -  3/4 
C  -  1 


0.778 

1.750 

3.263 

6.185 

414.36 


4.096 

8.858 


113.9 
9.22  x  10* 
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ANALYTICAL  EXPRESSION 


Laplace 


f„(s)  -  7p 


!  lx| 


Welbull 


Special  Cases 
of 

Welbull 


fxOO  -  f  (x  -  Y)6  “  1  e 


a  -  tra  +  |)  -  r2a  +  |)i  2 


y  -  -a  B  r(l  +  i) 


6  -  0.5 
6-1 
6  -  2 
B  -  3 

g-4 
6-5 
6-6 
B  -  7 
6-8 
6  -  9 

6  -  10 


u(x  -  Y> 


-0.09 

-0.26 

-0.39 

-0.48 

-0.56 

-0.63 


87.72 
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PDF 


Gamma 


ANALYTICAL  EXPRESSION  °3 


£x(x) 


(X  -  c ) 


T(a  +  l)b 


a  +  I 


e 


x  -  c 

b  u(x  -  c) 


_  JL 

b  -  (a  +  1)  2 


2 

/a  +  r 


c 


-(a  +  1) 


1/2 


3Ja  + 
a  +  1 


Special  Cases 
of 

Gamma 


A 

m 

1 

1.41 

6 

a 

m 

2 

1.16 

5 

a 

m 

3 

1 

4 .  S 

a 

m 

4 

0.89 

4.2 

a 

m 

5 

0.82 

4 

a 

m 

6 

0.76 

3.86 

a 

m 

7 

0.71 

3.75 

a 

m 

8 

0.67 

3.67 

a 

m 

9 

0.63 

3.6 

a 

m 

10 

0.60 

3.55 

-V  _ 


APPENDIX  D 


ANALYTICAL  EXPRESSIONS  FOR  STANDARDIZED  CDF'a 

CDF  ANALYTICAL  EXPRESSION 

o  ,  *  <  -vT 

Onifora  Fv(x)  -  —  +  i  ,  -/T"  <  x  <  /T 

X  2/y  2  -  - 

1  ,  /T<  x 


0 

• 

x  <  -/T 

+  jl  + 1 

12  /6  2 

» 

-i/K*  ^  x  0 

Triangular 

Fx(x)  - 

2  . 

-  2L.+  JL  +  I 

12  /T  2 

9 

o  <  X  <  /5* 

1 

9 

/T<  x 

T 

Gaussian 

V*>  • 

erf  (x)  •  ~  [*  e 

m  i. 

2 

dc 

Laplace  Fx(x)  " 


9/ 


CDF 


ANALYTICAL  EXPRESSION 


■  ..cial  Cai,e3 
cf 

Ganuia 


Veibuil 


a  “  0;  F  (x)  *» 

A 


a  *  1;  f’x(x) 


a  *  2;  Fx(x)  = 


a  *  3;  Fv (x)  = 


a  -  A;  f  <x)  * 


r 

x  +  1 

(1) 

r<i) 

r/Tx  +  2 

(2) 

r(?.) 

r/y 

X  +  3 

(3) 

ro) 

r2x  +  4 

(4) 

r(4> 

r/7 

x  +  5 

(5) 

u(x  +  1)  -  II  -  e~^X  +  ])}u(x  4  l 


u(x  +  /?) 


u(x  +  /’S') 


‘■(S) 


u(x  +  v?) 


Fx(x)  **  (1  -  e 


_  (x  -  Y) 
a 


B 


]u(x  -  Y) 


a  -  ire  1  +  |)  -  r2(i  +  |)]"8/2 


Y  -  -a1/6r(l  +  i) 


_  (x  +  0.447)1/2 

'8  -  1/2;  Fx(x)  -  [1  -  e  0-473  ]u(x  +  0.447) 


V  > 


APPENDIX  E 


TABLE  FOR  KOLMOGOROV -SMI HNOV  t.-<T  19] 


Signif icance  level  (a) 


.20 

.15 

.10 

.05 

_ jn__ 

900 

.925 

.550 

.975 

.995 

684 

.726 

.776 

.842 

.  929 

565 

.597 

.6<»2 

.?0a 

.  £  2 '  ’ 

494 

.525 

.564 

.024 

.  734 

446 

.474 

.510 

.563 

.66? 

410 

.436 

•  0 

.521 

.  616 

381 

.405 

.438 

,486 

.577 

358 

.381 

.411 

.457 

.542 

339 

.360 

.3S8 

.432 

.  514 

,322 

.342 

.368 

.409 

.48c 

,307 

.326 

.352 

.391 

.468 

295 

-313 

.336 

.375 

.450 

,284 

.302 

.325 

.  361 

.433 

,274 

.292 

.314 

.349 

.418 

.266 

.283 

.304 

.233 

.404 

,258 

.274 

.295 

.228 

.391 

,250 

.266 

.286 

.318 

.380 

.244 

.259 

.276 

.309 

.270 

,237 

.252 

.272 

.301 

.361 

.231 

.246 

.264 

.294 

.352 

.21 

.22 

.24 

.264 

.32 

,19 

.20 

•  22 

.242 

.29 

DECILE  BOUNDARIES  FOR  SOME  STANDARDIZED  CDF's 


Degrees  I  Significance  level  (a) 
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